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ABSTRACT 


To study criticality condition of Natural reactors of 
Oklo a computational model, ■which, describes neutron balance in 
ore, is developed using thermal neutron diffusion equation for 
homogenious medium. Criticality equation is modified to includ 
contribution from epithe^maal fission of U--255. Also included 
are teoaperature dependence of various quantities. Main parameti 
determining criticality are discussed in the background of geo- 
chemical information available about reactors. It is observed 
that moderation feature of Oklo rocks were as good as heavy wate 
and neutron leakage was much less compar-ed to man-made reactors. 

To study sustenance of chain reaction once criticality 
was achieved, bur nup /buildup equations of tvjolve important 

isotopes is solved by numerical integration scheEie. In most 

1 b 

favourable situations reaction would have stopned after 10 n/ci 
19 2 

to 10 n/cm flueiico. To sustain it further, hypothesis of 

continued mineralisation was examined for various criticality 

situations. Quantity of additional mineralisr.tion needed is 

21 2 

obtained for fluenco to reach 10 n/cm It was concluded 

20 2 

that either average fluence reached was 5x10 n/cm or if 
it was not so hypothesis of accretion of urrujium during reactor 
operation alone can not explain sustenance of reaction. 
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INTRODUCTION 


Once upon a tine seme 1.8 11111011 years lack there 
worked atleast six nuclear reactors underneath one quite 
corner of earth which today Is called ly the name of OinLo. 

They worked at the time when earth •■’as inhabited not ly men not 
ly fishes nor evr-.n ty plants but ly numerous tiny creatures, 

blta - 

virus and bacteria. One of the bacterial species called l«3re- 
green algee had ly then turned out Victorians in its promissii 
struggle for transforming whale of earth making its atmosphere 
oxygenic so that in future plants, animal and men could inhali 
tate earth. It was thoir labor that had made possible spring! 
up of these reactors, though it was unintended (see page 15)* 

Not before June, 1972 had man known about these little 
wonders. P.K. Kuroda says " Prior to this man believed that 
world's first nuclear chain reaction was sirtificial and had 
occur ed at the University of Chicago on December 2, 1942. The 
discovery of the OleLo phenomenon thus contradicted this view, 
since it contained experimental evidence that nature, not 
humans, had produced the first chain reaction. "[l] . This is 
from a man who predicted such occurancoo 16 years prior to 
the discovery of OldLo [2,3]. The story of discovery itself 
is a unique one. On June, 1972, a team of the French Atomic 
Energy Commission, Banzigue et al. [4], were controlling on a 
routine basis the input of uranium at the I'rench urani'um gas 
diffusion Plant of pierr elate when a sample of secondary 
urani\3m showed U-'235 content of 0.7171 ± 0,001)( , which is 
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sli^tly different from its thou^t to De constant value of 
0.7202 + 0.0006X on earth. At first it was thou^t that 
this minor deviation is because of experimental error but 
repeated experiments led thCT to rule out this possibility 
and suspect contamination by commercial depleted uranium in 
supplied uranium. Search for the source of contamination 
began, it had been’ supplied by French Comp'any, CQ'dUF, opera- 
ting two mines of OlfLo and Mounana in Gabon. Soon it was 
known that uranium traded by this company for one and a half 
year, amounting to 700 tons, had less than normal U-255 
concentration. On examination of minerfiL of Oklo to their 
surprise ore itself was found depleted in U-235. In one rock 
sample depletion was up to 0.4417- • Various hypothesis were 
advanced to account for this anamoly. If it was because of 
nuclear chain reaction their would be abundance of fission 
products. Soon ro.ro-earth fission products (Nd) in large 
quantities were detected by mass spectremetry . This was a 
decisive evidence for nuclear chain reaction and in September, 
1972, Feuilly et al. [5] confidently announced that Oklo 
deposits had been the site of Natural chain reaction. By 
larch, 1973 CEA formalised investigation of Oklo phenomenon 
as " project Franccvelle” under leadership of Dr, R. Neudet [6] 
thus initiating on site examination of mine. In the beginning 
:his team identified three distinct uranium rich zones in 
•oefewhere marked depletion was found. One record sample 
iontained surprisingly low value of 0,296/. II-235* One sample 
iven showed slight enrichment with U-235 content of 0.74/* . 



Fortunately good correlation between U-“23‘5 isotopic abundance, 
uranium abundance end fission products abundance was establish^ 
indicating that reactor fossile had remained intact after 
reaction stopped in it. Rocks of Oklo had kept complete record 
of a unique event. Continued investigation soon led to identi- 
fica-tion of six fossilised reaction ^ones in same site. At tha 
time in a reflective mood Neudet wrote ’'we don't know which 
should be considered the more astonishing, that these (j/eactor) 
mechanisms existed or that we find such detailed evidence 
nearly two hillion 3, oars later" [6]. Geological age of working 
of these reactors by that time had been oetimatod to be closo 
to 2 billion years. 

Since then two international meetings on Oklo phenomenon 
have boon held with patronage of IAEA at Libreville in June, 

1975 [7] and in Po,ris in September, 1977 [8], Wide range of 
extensive multidisciplinary studies were reported in these 
conferences ranging from geological, geociianical, radiological, 
patrological to i-c-actor physics aspect of Olelo phonomenon. 

These studies have illustrated how Oklo phenomenon challenges 
over understandiiig of wide variety of complex problems like 
ageing of geological and gc-ochemical froaturcs over a long 
time scale, stability of radiation dimnagi, features in insulators 
radiation induced corrosion, ccmpositlon and circulation of 
deep ground water, redistribution and stability of elements 
in nature over a „vory long time scale. Also Oklo phenomenon 
is an challenge to analytic methods used in variety of fields 
like earth sciences, irradiation effects in solids, radio^ 
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chemistry and so forth. Only the most refined methods used 
in geochronology, rare-gas determinations and ion-microprobe 
analysis can be exploited in the "non-closed" geochemical 
syston of the Oklo reactors, where in addition to other processe 
neutron fluences have been active to further modify isotopic 
abundances. Besides these probligns reactor physics aspect of 
these reactors is challenging. Problaia of understanding 
sustenance of reaction, operating mod e of reactors and duration 
of reaction are still open. 

One of the discovery of far rcachirig consequences from 
Oklo has been a finding that quartz variety found in Oklo 
ore has the hipest track retentivity ever reported for an 
insulator [9]. Laboratory experiments sb ov;- that fission track 
registered in the Oklo quartz should be stable at 300 ^0, over 
a time scale of 10 billion years, and this is indeed a world 
record for solid-state track detectors. Freshly induced as 
well as fossile tracks are still stable at a temperature of 
1300^0 vrhich is only slightly lower than the fusion temperature, 
1500 °C, of Oklo quartz. Very pure varieti es of synthetic 
quartz grown by tho fusion zone technique resists annealing 
of radiation damage effects only till 40C°C. Researchers had 
to postulate that much higher stability of Oklo quartz is 
related to peculiar’ impurities that act as stabilizing agents. 
This is a fine example of how man in his quest for bettorhood 
was loyal to purity and missed something, its quite natural as 
ideas are pure and contain deniel of plurality. 
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Another interesting study is by A. I« Shlyakhter [10], 

who developed direct test of the constancy of fundamental 

nuclear constants from measured isotopic composition of Sm and i 

in samples from Ohio rea.ctor core. Prom fission yields of Sm 

and Eu isotopes which had got irradiated by an independently . 

determined fluence of reactors he calculated their capture 

cross-section values some 1.8 billion yci'.rs ago. These isotopee 

were chosen because of there low lying neaj:‘ thermal resonidcos 

which would shift if strong intera.ction coxistant, g , or weak 

s 

interaction const-rao t, g , or Couluinb coupling constant, a, 

w 

were to change with geological time, thus shifting capture 
cross-soction vaiues. From calculated value of capture cross- 
section considering error involved in their determination he 


showed these constants 

are constants within limits given bolow 

1/ gg 

dgg/dt 

1 

H 


1 /a 

da/dt 

(yr-^) 

10-” 

1/Sw 

dg^/dt 

1 

H 

; 2x10 


American’s have shown their interest in yet another 
direction [11, 12, 13»14]. In their assessment, Oklo phenomenoi 
is an excellent exaiuple of confinement of fission products, more 
popularly called radioactive waste these duys, which had not 
migrated out of reaction site for around period of 1.8 billion 
years. Except fisoiogenic rare-gasos practically everything 
has remained confined in Oklo reactors for 1.8 hillion years. 
This they think has r. direct bearing on tho pressing problem 
of long term disposal of man-made radioactive waste. 


Oklo for 
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them testifies safe storage of radioactive waste in natural 
geological setting is feasible. Study of geochonical features 
of matrix in which fission products remained intact in Ohio 
could provide insight so that ways can be explored in which 
such features can be inbedded in design of nuclear waste 
storage facilities to minimise hazard. 

As far as reactor physics aspects are concerned one 
extensive review [15] and six reports [16, 17, 18, 19, 20, 21] 
have been published in above mentioned tv/o international 
meetings. These psipers have dealt with criticality conditions 
and problem of control and sustciinance of reaction. Major 
inclination in them have been favourable to the hypothesis 
of sustainance through burning up of initially present poisons. 
Our work banks on the hypothesis of reaction sustainance 
through continued mineralisation. ¥e have undertaken to 
examine consequences and ground of this hypothesis. 

In first chanter review of geological, geochemical and 
physical feature^ of natural reactors is done. Also described 
is the formalism adopted to represent physical features of 
reactors. In second chapter is developed neutronic model for 
reactors. Conditions under which Oklo reactors would have bee a 
critical are discussed. In third chapter isotope burnup /buildup 
model is developed to study sustenance of reaction. Quantity 
of additional accretion of uranium is computed for sustenance 
of reaction till desired fluence. In cliapter four is discussed 
general possibility of occurance of chain reaction in nature. 
Also discussed is canclusion and limitation of this study. 


Chapter - I 


GEOLOGICAL, SEDIMENTOLOGICAL AND PHYSICAL FEATUdEG OF OKLO REACTQPS 
1^1 On possibility of nuclear reaction in uro.niitm ore: 

In 1956 Zuroda[;^,3] for the first tine visualised 
possibility of naturaJ. formation of critical assembly in mines 
of uranium. Analysing, J ohanngeogenstadt nitcbblendo he computed 
e, p, p and f as a function of geological time and established 
that for water to uranium ratio greater than half uranium mine 
would be supercritical before 2100 million years ago. Through 
his study he thought he had explained why ages of large scale 
uranium deposits never exceed 2100 million years as any such 
deposit at that time would destroy itself. 

Discovery of occurance of large scale fission in Oklo 
uranium mine confirmed his visualisation btit falsified his 
conclusion, as the fossil e of reactor site is still intact and 
there is no evidence of destruction. And here precisely lies 
the mystery that how i/as nuclear reaction at Oklo sustained and 
controlled once criticality was achieved. As to how criticality 
was possible is in I’ough terms evident from Pig. 0, which gives 
content of uranium -235 in natural uranium as a function of 
geological time. Relative content of uranium-235 was higher 
in the past as decay rate of uranixmi-235 is higher than that of 
uranium- 238 . 1465 million years ago uranium-235 enrichment 
was 2.3/. . This would mean Tarapur Reactor could be fueled by 
natural uranium at that time. Also any assemblage of concen- 
trated uranium could go critical with natural water around or 
before that time. 



99 **/ 
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To go into details of working of enearthed natural 
reactors we will have to trace origin of Okie deposits and 
subsequent evolution leading to their working. As the physical 
history of reactors ctm not be separated from geological history 
it would require review of geological and geochemical aspects 
which were central to unfolding of this phenomena. 

lo2 Geological and geochemical features : 

In Table I results of several studies regarding age of 
uranium deposition have been ccanpiled. Average of these esti- 
mations comes to 1905 + 151 million years. The reactors must 
have worked during or after this time. 

Oklo mineralisation geologically was an terrain phenomena 
a thorium content in uraninite crystal is very low. This is 
only possible if deposition took place from w ery low tonperature 
(around 70°C was surface temperature at that time) aqueous- 
solution because tlioriim is insoluble in low temperature aqueous 
solution [22]. Also preservation of fossil algae in the host 
Irancevelle series is indicative of that [23]. 

About 2000 million years ago land mass over earth looked 
quite different. South America was close to Africa and Oklo 
was situated at the ridge of two evolvj.ng continents. Fig. 1 
shows position of Oklo in the global map at the time of ore 
deposition. 
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TABLE 1 


Age of Mineralisation 

: (t;) 

' in Million years '<• 

. , , .. , < 1. , Tl I T , 1 

Method of study 

1 

f 

: 1800-2000 

1 

K/Ar isotopic study [42] | 

' 1750 

U/Pb isotopic study [43] 1 

; 1800-2050 

U/Pb isotopic study [44] i 

, 2000 

1 

]ld,Sm,RUj,Th isotopic [44] i 

1 

abundance study ! 

1600-2100 

! 

Parameter fit model [41] ; 

of Nd isotopic abundance ' 

1800-2100 

1 

For different Thermal [45] ' 

i 

fluences published j 

145 ! 

Ed isotopic study ; 

. 2000+100 

From Uranium fission [46] ; 


d ata i 

, 2050+30 

U/Pb study [47] * 

. 1730 

For motion age of [48] | 

quartz crystal from 

fission track analysis j 

i 1700-1800 

( 

1 

K/Ar and Rb/3r isotopic [49] ; 

analysis j 

1800 

K/Ar and Rb/Gr isotopic [42] ; 

analysis | 

; 1930 

Th/tl/t>b isotopic analysis[51] ! 

^ 1900 

1 

! 

; 1980+500 

1 

Ed isotopic study [51] j 
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African continent is made up of four consolidated and 
granitized Cratons. The West African, Con^.olese, Kalahari and 
Nilotic Cratons are represented as region A, B, C and D respec- 
tively in Fig. 2(a). With these Cratons are associated internal 
"intracratonic" basins with shallow detritic sedimentations and 
int ercratonic zones with their own deep sedimentations. Intra- 
cratonic zones are located at the periphery of the cratons on 
the fold axes and are called 'mobile zones'. Gabon lies in the 
north-west part of the Congolese Craton. Franc ev ell e basin 
which is in south-east of Gabon, and in which lies Ohio is one 
of the intracratonic basins of the Congolese Craton and is marked 
22 in Fig. 2(b) . 

The present structure of Africa only got established by 
the end of precambrian era (500-600 million years ago) [24]. 
Before that geological history of Africa is marked by series of 
brief paroxysmal eveuts in prolonged periods of stability. 
Surfacial geological cj/cle of erosion, transport, sedimentation 
and folding is called Orogeny. The final paroxysmal phase is 
called orogeneses. Chronologically in Africa there have been 
four major orogeneses ; Precambrian A (500-600 to 900-1200 
million years ago), Precambrian B (900-1200 to 1800-2000 million 
years ago) , Precambrian C (1800-2000 to 2500-million years ago) , 
Precambrian D (before 2500 million years ago). The Franc ev ell e 
basin can be assigned to the precamberian B orogeneses as age of 
deposition for it falls in that period. 




fht. 2, rc!-olo5^Jc-?l Afrirr;. <Jratons. A » West African crafon, 8 « Congolese 

crc/^nn, C- ::r.'!tcn. S P" riripr-si r/ 

by nurr>tcr?rs. (fhjru T: Mohn<. cr-:J J. C B0$'0f^S, Tim (Jklo 
£nergy /•ty>^r:-^-'0s4 .^975) 



14 


Through an independent study also one can extablish that 

uranium ore formation in sediments could not be older than 

precambrian B orogeneses. Before 2200 to 2000 million years ago 

earth's atmosphere had little oxygen. Because of lack of oxygen, 

surfacial uranium in volcanic ashes and Basaltic rocks was hi^ly 

44 “ 

immobile as uranium existed in U chemical state, in which it 

^ 4 “ 

is insoluble in water, and only in its oxidised state U it 
is soluble in water. Brom chemical immobility of uranium in U^'*" 
state we can infer that uranium existed in very lovr concentration 
(about 10*"^/. by weight) in volcanic ashes and older Basaltic 
rocks. Around 2200 to 2000 million years ago rise of new form 
of becterial life called blue-green algee, which was the first 
organism capable of photosynthesis, liberated oxygen in large 
quantities. Thus making earth's atmosphere oxygenic also 
increasing content of dessolved oxygen i3& >,ater. This was not 
only important for future evolution of life but also played 
irreplaceble role in oxidising uranium to state enhanced 
by the process of surface leaching, by the process of surface 
weathering because of torrential rains, and by rock erossion 
because of winds. This only made possible mobility of uranium 
as Uranyl ion, UO 2 , along with metals having similar geocheaical 
behavior. With its host, water, it move around surface of 
earth along rivers, infiltrating sediments till it got reduced 
in geological traps back to U^"*” state [9 ]. For uranium geolo- 
gical traps in sediments were locations with hj.^ reduction 
potential where uranyl ion got detained from moving water. If 
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conditions had been such that water flowed through traps for 
long time uranium concentration in them would increase to make 
those sediments uranium ores of future. These series of phenomena 
not only explain formation of uranium minerals but also explain 
why there is no large scale known surfacial deposit on earth 
older than 2100 million vears. 

Existance of concentrated localised carbonacious matter 
in sedements could have been ideal reduction medium for uranium. 
And fossilised life colony or organic ooz in newly formed 
sediments cross-cutting waterbed probably provided such condition 
for formation of uranium ore in Prancevelle basin. Pour such 
large locations illustrated in Fig. 3. have been unearthed in 
Prancevelle basin, those of Boyindzi, Mounana, Oklo and Okelo- 
bonde. All of them are situated in a same flexural zone effected 
by sliding. First two, at the contact of basement are subject 
to break deformation, while the other two are implicated in 
extensive flexible structure resulting from s-ame sliding 
tectonics. Oklo uranium deposit is situated 1500 m south of 
Mounana deposit. Concentration of uranium is these large zones 
would 'nt have gone more than 0.4/. to 0.6/. by wei^t in primary 
mineralisation [25]. Though there would have been little 
variation in degree of mineralisation due to locally variable 
sedimentological conditions. Geological studies indicate that 
primary mineralisation stopped because of continued sadimentation 
as mineralisation zones reached depths inaccessible for uranyl 
rich Tiater to penetrate it [26]. 






or: » 


UrrJilu^fi 


ri in feral 1 £?.:.■' ■■/iCiii' in B'rfiXiceveTle Basin 
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Upper figure in Pig. 4 maps zone of primary minerali- 
sation in Oklo projected on horizontal plane as seen vertically 
where as lower figure is a sectional view along vertical planes 
raarked A-B, in upper figure. In Franc ev ell e tasin there has 
been two superposed sandstone-pelite cycles. In the fore deep 
basin uranium is localised in a 4 to 10 m thick layer in upper 
part of basal sandstone. Maximum size of mineralisation zone 
is 900 to 600 m. Figure 4 illustrates situation of deposit 
today. Actually uranium mineral went through long history of 
tectonics and further geochemical rearrangements before aquiring 
this shape. One u'jique feature of Oklo as chi-stinct from other 
three deposits in Francevelle is that it is close to basal rise 
which conditioned tectonics and influenced sedimentations. Fig. 5 
illustrates how priiiiary flat mineralisation zone because of 
future tectonics under?,; ent change. This ufhole process of pre- 
Cambrian B orogenesis must have taken millions of years before 
leading to present rock shape. These tectonic movonents lead 
to formation of extensive shear trouts in Oklo mineralisation. 
Whereas in Boyindzi and Maunana sedimentary rock got fractured 
and in Okelebonde little sliding occur ed. The formation of 
shear troughs played unique role in the formation of local 
t ectono-sedimentary traps in primary mineralisation zone which 
lead to reconcentration of uranium in those areas after mineral 
sediments reached deep burial depth. Pressure and temperature 
conditions at that depth coupled with paroxj^'smal tectonic activity 
lead to total desilicif ication of sandstone by oxidising currents 








5- folding of rxineralisation zone by tectonic activity 
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in shear troughs and uranium got concentrated in extension 
zones forming anamolous uranium rich lodes within large normal 
mineral zone. When this process of concentration of uranium 
was going on at saae point eventually criticality was attained. 
This silicon diagenesis lead to compaction and settling of 
earth to the degree in which we find Oklo ore today. Evidence for 
this diagenesis is substantiated by patrographic studies [27] 
which indicate two distinct type of uranium ore in Oklo x^ith 
distinct modes of formation. In normal ore uranium appears 
mostly as micr o-veinl ets of pitchblende, with corresponding mean 
uranium concentration about 0.5/. , embedded in a matrix of clay 
mineral rich in organic matter. In deposition of this typo 
of ore pred eminent rolo was played by organic material. In the 
anamolous uranium rich zones sandstone is predominantly chlorite 
rather than clay and uranium appears as a cluster of tiny 
( 100 Mm) and x^^ell crystalised uraninite grains and uranium 

shows about 50 fold increase in concentration over 0.57. value 
in normal ore. Certain spots even have uranium concentration 
upto 60/. . Migration of oxidising solution could only be the 
main source of reconstruction. And this would need moving out 
of silica to make place for uranium. 

Six sites of these compact accumulation of very rich 
uranium ore has been identified in Oklo, All of them virtually 
lie in a single block flattened in the direction of stratum. 

Average thickness of these lodes is about one meter and their 
length about 10 to 20 meters. i'ig.s.'’6 shows horizonto,! 
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projection of these zones maped according to the radioactivity 
counts in radiation counter, indicating uranium concentration 
today. Fig. 7 is a sectional view of these compact zones along 
vertical plane A-B, C-D and E-F. Total surface area of these 
zones is about 1.5/ of surface area of primary mineralisation. 
Average uranium concentration of these zones is greater thcai lOXby 
wei^t today. It iras only in these zones that marked depletion 
of uranium-235 is found. This would imply that nuclear reaction 
took place onl3/ in than. Certain tectonic episode coupled with 
continued increase in bur i el depth had triggered their formation 
within primary mineralisation. This was made possible because 
of uniqueness of rock structure of Oklo deposits. Fig. b illus- 
trates fold structure of Oklo which stroches 500 meter in ITorth- 
South direction and reaction zones are shoun at base of these 
folds lying in the extension zone of Dome slightly away from 
the main slippage. Rise of this Dane was central for formation 
of uranium rich lodes and attainment of cri'fcicality in them. 

When reaction began it resulted in local thermal motamorphism. 
Thermal syphoning progressively expelled water and organic 
matter from core of the reaction zones. In this process coarsely 
crystalline uraninite was formed and eventually burning up of 
uranium brought reaction to end. After finishing up of criticality 
settlement of earth took place once again in which carbon played 
dominant role that is why this phase is called carbonate diagenesis. 

Geological history since primary uiinoralisation has been 
illustrated in Pig. 5. Depth which mineral sediment reached has 
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Fig„ 9- Geological history since primary mineralisation 
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been obtained from working ont of tectonic history and also by 
fission track studies [9]. It was in the later part of desi- 
licanisation process that urani-um got concentrated enou^ to 
start chain reaction. At that stage reactor wa.s buried 4.8 Ions 
deep. Generalised T,,riation of lithosta,tic and hydrostatic 
pressure with depth has been plotted in fig. 10. Dotted line 
indicates pressure condition at the finish of d esilicanisation 
process. In Fig. 11 is plotted pressure’ temperature geothermal 
gra,dient at Oklo. lino 3 is an isochore oi the fluid in rock 
contemperary to the formation of primary vicn of quartz at tho 
begining of desilica diagenesis. Line 2 is an isochore of fluid 
contemperary to tho formation of final quartz at the finish of 
nuclear reaction. Line 1 is isochore contemperEiry to tho 
carbonate diagenesis after reaction site cooled down geologically. 
From these curves one can infer that temperature of fluid at 
the time of establishment of criticality must be between 180 °C 
to 230°C. If surface temperfiture in prcccinbrian times was more 
than 50*^0 as one study indicates surface temperature to be 70 °C 
to 80^0 than temperature of sediments would be 210*^0 to 250 °C. 
Gross-sections used for studying neutronics have to be modified 
for this temperature. 

I.J Uranium content in reaction zones • 

In Table 2 are comp ill ed content of uranium in reactor 
zones as well as amount of uranium- 23 5 depletion observed [35] in 
reactor sites today. Prom these figures simount of uranium-235 
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burnt if reaction tool place 1900 million years ago have been 

deduced. Also total amount and volume of materials in these 

zones have been tabulated with the assumetion that 30 % of 

material in reactor zones was uranium. Idealising reactor zone 

shapes by one meter thick slab with square surface, geometric 

Buckling has been calculated. Bor all the reactor zones it is 

sli^tly more than 10 m”^. Bor reactor systems 5-4 and 5-6 

the surface area will be much more than iderlised geometry area. 

This would imply geometric buckling would be higher than calculated. 

If reaction had began originally in smaller volume and grown to 

the volume observed today than again geometric buckling would 

“2 

be hi^er first and then will reduce to 10 m by closing of 
reaction. 

Measurement of content of uranium, U*/. , and partial 
density, d^, of uranium were carried out for various samples 
from ore site. Partial density is defined as weight of uranium 
per unit volume of sodementary rock. Result of the survey [28] is 
plotted in Fig. 15. I or analytic purpose best fit curve can be 
idealised by following equations 



= 0.52 - 0.48 X U% /lOO 


( 1 ) 


From above equation one can find out number density of 
uranium for any value of content of uranium ranging between 10/. 
to 50/. by weight. 




Partial Density in gms/cc. : 

Fig. 13 Relation of partial density of 

Uranium and content of Uranium by 
weight in Oklo rock 
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1^1 Content of rest of material (gangua) in reaction :iones : 

From the same survey content of rest of material in 
reaction site as a tunction of uranium content in rock and 
geological time is infered as [16]-; 

X = X exp (-0.17 5 t) (2) 

where X is defined as weight of gaxgue per unit weight of 
UO 2 in sane volume, 

and X is geological time in billion jofvs. 

Exponential terra comes from taking into account of 
transformation of uranium into lead by natural decay. Other 
terms come beoause gangue does not include water which is 
chemically bound with mineral also a factor comes because 
wei^t of UO^ is different from wei^t of uranium. 

Io5 Content of dijfer-.'iit constituents oi gangue : 

Silicon and Aluminium are major constituents of Oklo 
rock as clay and sand were host sediments where mineralisation 
took place. Magnesi-um, Iron, Manganese, Potassium and Calcium 
are easily traped by organic rich lences along with uranium 
from low temperature aqueous solution. These are the elements 
which constitute most of the gangue thougti there are thirty 
odd elements found in rock* Average of measuraaient of relative 
content of different elements in reactor zone. ? are as follows: 
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SiO^ 

42/. ; 

Al^O^ 

28% ; 

PeO 

10.5%. ; 

Fe20^ 

4,5% 

MgO 

7.5% i 

K^O 

2.5% ; 

CaO 

2.5% i 

Na20 

0.6% 

Ti02 

0.7% ; 

MnO 

0.5% ; 

BaO 

0.2% ; 

Others 

0.5% 


Prom neutrorics point of view content of important 
trace elements was found out to be following 


¥ anadium 

200-800 

P » 

Or, ITi, Ca 

50-100 

ppm 

Boron 

20-250 

ppm ; 

Sm 

10-30 

ppm 

Gd 

7-20 

ppm ; 

Li 

100 

ppm 

Hf 

20 

ppm 





For simplicity and to develop furthei* analytic scheme 
constituents of gangue can be divided by infei’ence from actual 
survey as given below [16]; 


Si02 

0 

« 

44.5 

- 0.05 

X U% 

AI2O3 

0 

50.5 

- 0.05 

X U% 

MgO 

0 

11.7 

- 0,08 

X U% 

PeO 

0 

0 

10.0 

+ 0.15 

X U% 

K2O 

0 

0 

2.5 



MnO 

0 

0 

0 , 0 + 

0.01 X U% 


These relation have been derived by research team headed by 
Neudet [16 ] from actual measurement of content of different 
element oxides in roch samples with different uranium concen- 
tration. Also these relations have been checked by study of 
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geochemistry of desilicanisation process. In above given set 
of equations contents of Si02> Al20^ WgO are a bit over- 
estimated to normalise but because of their low scattering c 
cross-sections these ov er estimations account for other elements 
neglected. MnO is olso overestimated but because of its high- 
capture cross-section it takes care of other elements neglected. 

F eO estimate is low as and F eO have been bracketed and 

their combined scattering can be better estimated by taking 
sli^tly lowFeO content. 

Number densities of constituents in gangue can now be , 

written in simple consistant analytic iOrm as function of 
content of uraniuifi and geological time. 

The trace clement in gangue which have hi^ neutron : 

capture cross-section have to be separately'' taken into account. j 
Macroscopic absorption cross-section of ali. of them can be 
represented as equivalent to sane quanti'by of Isotopic 

geochejmical stxidios indicate their ostimateci content in reaction 
zone could have been equivalent to 50 to 250 atoms per 
million atoms of uranium [16]. 

2 06 Wrxter content in reaction zone •" 

There is about 9% to 11/. water by weight in gangue 
which is chemically bound with gangue. Thore is also possibility 
of existence of water in the pores of rock. In fact hypothesising 
existance of water in rock per os when sediments 'were buried deep 
is vary essential to the thesis of desilicanisation and decorbonat< 
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diagenesis- During primary mineralisation, porosity planes, 
where mineralisation took place, must h£-ve been saturated 
with water. As the buriel depth of sediments increased there 
was possibility of water boiling out as precambrian surface 
temperature was about 70 ^ 0 . In that case whatever water remained 
in rock pores could htvo a.t deep buriel dcpiiis moved into shear 
troughs formed by tectonic activities and it was there that 
reconcontration of uranium took place. It seems quite acceptable 
that rock pores in reaction zones must have been filled with 
water. Estimation of volume of rock pores is difficult. But 
generally 6 to 12 '/. volume of rocks is microfractures for 
old sedimentary rocks. Two kind of water in Oklo rock can be 
formalised in following way to express content of water. 

Y = w eight of water/volume ^ y , 

" w eight of trb 2 /'^ olum e 

whore m* is acontent of water chemically bound to gangue 
per unit weight of gangue. Xm ' gives weight of chemically 
bound water per unit weight of UO2 foi- unit volume, m is 
weight of water in pores per unit weight of UO2 in unit volume. 

Wei^t of loose water in 100 gns of rock is 
= 1.34454 X U/. X m 

If d was density of water in those condition and 
w 

d^ was density of rock then volume percentage of loose water 
is 

d^ X m X U’/. 


V/. 


1.345 
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This has to he always less than volume fraction of 
rock pores. Wow expression for water content becomes 

d xV /C 

y = Xm' + 0.744 where [6/. to 12'/] 

Wate^r is pores during and after dc-silicanisation process 
was subjected to pressure between 800 bax-s to 1000 bars. and 
its temperature was 150°C to 300^0. Water in this region is 
supercritical and its density varies smoothly with temperati’.r e[40 
as shown in Tig. (l-l). In the region of interest straight line 
approximation can bo made to analytically formulate tsnperature 
dependence of loose xvatcr density. 

d^^ = - 9.4 X 10"^ T„ + 1.106 

w m 

where T^^ is temperature of milieu in degree centigrade. 

Wow expression for T becomes 

Y = Xm‘ + 0.744 (-9.4x10“'^ + 1 .106)V/./(U'/. d^) 

(4) 

Here we have formulation of T wh: ch is sensitive to 
measurable and physically constrained parameters Y'/, , d^, m'. 
m' was averaged over reaction zone 2 and wus found to be close 
to 0.10 (lOX of gangue). Density of dry Oklo rock has been 
estimated to be bet'-een 2 to 5 gms/cc for varied uranium 
content. Fig. (15a,b,c,d) shows variation of Y with V*/., U'/ , 
m’, d^ at temperature 200°C 1.8 billion years ago. With most 
optimistic combination of different parameters Y coiOLd not 
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have been greater than 0.5. If we take volume fraction of 
loose water to be about 10/. (= V/ ) and ui’.xaium content to 
be around 24/. (= U/ ) as v/as in most of Cklo reactors, Y turns 
out to be less than 0.32 for = 2 and le>ss than 0.28 for 
d^ = 3. By simply averaging densities of components of rocl: 
wo can estimate that density of rock is cJoser to 3 gms/cc 
and we can say with seme certitude that Y could not have been 
greater than 0.3 in Oklo reactors. For higher \jranium contents 
it is evon less. 



CHAPTER - II 


STUDY GF CONDITIONS OF CRITICALITY 
2ol Reactor physics model : 


In this chapter reactor physics model for hatural 
reactors of Oklo has been constructed vrhich will enable us 
to assess conditions under which natural asseniblage of uranium 
could go critical. I' or evaluation of criticality conditions 
we would work with the assumption of homogenised medium 
disregarding spacial distribution of uranium and other 
materials in Reactor zones as well as whatever lumping of 
uranium exists in them. Also we would work with the assumption 
that all material in reactor zone is exp.osed to same neutron 
flux. 


For any assemblage of fissile material to sustain 
neutron chain reaction effective neutron mi'iltiplication factor 
(Keff) should be one or greater than one. Prom hcmogenious 
thermal diffusion equation and Permi Age theory expression for 
Keff is given by foilovring equation [36]. 


K. 


K 


ini 


ef f 


(1 + l^B^) 


exp 


( 5 ) 


where 

is gecmetric buckling of assemblage 
is diffusion area for thermal neutrons 
is Age fron fast to thermal energies 
Kinf is infinite multiplication factor. It is an intrinsic 
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property of materials in Reactor and is riven by the product 
of four factors. 


K 


inf 


p • f • e ‘p 


( 6 ) 


Four factors 
p is 
f is 
£ is 
and p is 


are: 

regenerration factor 
thermal utilisation factor 
fast fission factor 
resonance escape proi ability . 


Now we will elaborate on different facton-s constituting 
Eqs. (5) which chai'ac terise behavior of neutrons in Oklo 
reactors. 


2^1„l Regeneration factor, p ; 


It is defined as number of neutrons produced per neutron 
absorbed in fuel. Definition of fuel is to some extent arbi- 
trary. Thoug^i in OIlIo uranium-235 was the only fissile isotope 
but we will define fuel as both uranium isotopes because both 
isotopes of uranium are chemically closely associated. Then 
regeneration factor is given by following equation. 


where 


p = 


N 


25 


N 


28 


4 


, il 28 28 s 
a 1. 25 a 


is number density 
is number density 


of uranium- 23 5 
of uranium- 23 8 


( 7 ) 
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25 


of urani'uia-“255 


25 

a is a'bsor;''<tion 

8- ®®~seotion of urani\Ain-235 


28 

a is absor-Dti on rv-v^^ 
a ^ « cross- 

NorA'^or varies with 


Section ox xiranium“238 


2® 25 S'eological time as decay rates 

of two isotopes are '■3.ifferen + 

Decay ecivations of two isotopes 


are 


N 


28 ^^28 ®^P^Q-153691 x) 


,0 


**25 = ®25 (0-972156 x) 


where N^o Npcr sxe nximhen -i.- j. -i 

28 25 densities of two isotopes today 

and T is aeologics,! time i'-o ■ . 

^ of billion years. It increa- 
ses as we go in past. Relation ^ •,.+ 

ab-undancc of two difiereiit 


isotopes is universally constant 
^28^^25 ~ 137.85. This impli^g 


at any given time. Today 


= 137.85 (-0.818465 x) 


(8) 


From above equation 
was more in past. Fig. (i) 
Substituting above equation 
geological time 


Can See that uranium- 235 abundance 
plotted fi'om above equation. 

(7) we get n as a function of 


2ol.2 Thermal utilisation factor , f 


Thermal utilisation - , _o 

factor IS defined as the ratio of 

mmiber of thermal neutrons absorb^^} -rv • j. a. - - +o 

'-'-vDea per unit time m fuel to 

I 



the total mm her of thermal neutrons absorbed per unit time 


in reactor. For ofiectively homogenious reactor it is given 
by following equation. 


f 


(1 + 


n 


i=l 



4 



(9) 


where 


a 


is ther'asl neutron absorption cross-section 
ijln 

of i Diaterial in reactor zone. 


N- is number density of i "■* material in reactor zone 


a 


is absorption cross-section of fuel = 


28 ’■"i 

^ 28 ^a 


and N. 


N 


28 


+ h 


25 


From Bqs. (1), (2), (5) and (4) number densities 
as a function of uranium content U/. , geological time x, 
neutron poison content x, and volume content of water in 
pores VXsCan be evaluated. Substituting tliem in Eq-(9) we 
will get f(U*/. 5 X, Xy V/.. 


2^1^^ Resonance escape probability, p ' 

This factor talces care of the absorption of neutrons 
in epithermal energy range- Uranium is main absorber of epi- 
thermal neutrons au.d other materials are weak absorbers. 
Expression of resonance escape probabilit;>' for every element 
is 



Pi 


exp(- 




f 2- / 

^ s/e 




B, 


■eff 


X 


i dE' 


E 


th 


1 + 


( 10 ) 


where 


is average logarithmic decraient 

n 




i=l 


f / N . . 
s/e 11 



1 


i=l 


s/e 


TvT 

1'* jj 


■'s/e 


is total macroscopic epithermal scattering 
cross-section 

= f, B. aS 

1 s/e 



is microscopic epithermal scattering cross-section, 
of i"^^^ clement 



^eff 


is microscopic epithermal ahsoi-ption cross-section 
*t*Vi 

of i element 

is effective resonance absorption integral. 


Bor weak absorbers in our reactor like uraniuin-255, Si> 11, 
Mg, Be, Mn, K and 0^ above equation for each reduces to 


Pi 


exp(- 



^inf ,dil. ^ 


( 11 ) 


where I. ^ is infinite dilution resontaice integral, 

inf. dll. 

Epithermal absorption by hydrogen can be negd ected. 
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For uranium-?J38 N20 neglected, 

because of its large quantity in reactor zone. Experimentally 
determined values 01 integral in Eq. (10) are available for 
different val.ues of 28 /^ ^ /e ^ S' oll'^''f'r4ng; imperical 

relation can be used for integral. 


leff^28) = exp ( 0.46 log( 1.123) 


V 


28 


exp ( ■ 


N 


28 


^eff^ ' 


( 12 ) 


5 ^s/e 

Using equations (11) and (12) total resonance es.cape probabilit; 


n 


r = Pja (,U, Pi ) 


( 13 ) 


p is a function of U'/ , t, x, V*/. 


2= 1.4 Fast fission factor, e 


It is a measure of neutrons born in 11-238 fission by 
fast neutrons. In liomogenious reactors it is quite close to 
one but in Oklo reactors fast effect is considerable because 
U- 238 /water ratio is quite hi^,. Calculation for this factor 
is difficult without lot of assumptions. ¥e assume that only 
elastic collision of fast neutrons with vrater brings it down 
below fast fission threshold of about 1 K eV and elastic 
collisions with other materials present in reactor do not alter 
fast fission probability by neutrons. Another assumption we 
make is that only inelastic collision with U-238 is possible 
and it brings down neutron below fission threshold. In another 



vrorda we have assumed for calculation of fast effect that Ohio 
reactors only have hcmogeniouo mixture of U-238 and water. In 
this case e is given by following equation. 

e = (1 + 0.69 N28/\)/(1 + 0.563 


where K is number density of water, 
w 

This would mean e is a function of x, U'/. and V/. . 

p 

2ol^5 Diffucion arca^ 1-2 : 


Diffusion, area is related to meaji square distance 
travelled by neuti'on in thermal energy range. It is given by 
follox'^ing equation 


^ ^th^ "^a/th 


(15) 


wh er e 


"a/th 


is total thermal macroscopic absorption cross- 
n 


ection = 


and 


D^n is tliermal diffusion coefficient 


th 


_ ^ - N / n 4 ^a/th . ^a/th i^^o 

l/(q 5 T:^ 


■)) 


1-p 


o 


is thermal macroscopic total cross-section 

= t ''i ''t 

i=l 

n is average of cosine of scattering an^e. 



Equation (15) will give as a function of UX? 
and Y’/.. In Oklo reacnior medium scatter irjj. cross-section will 
have some variatiOii in thermal range because of 
interaction between reutron and the molecular structure of 
medium or from interaction betvj'een neutron and the crystal 
structure of the mcodium. That is why calculated value of 1 
from Eq. (15) will be- an underestimate. 


2^1o6 Neutron age, t 


Age is related to mean square diGtance travelled by 
newly born fast neutron to thermal energy range. It is given 
by following expressions 




where 


u-u 


( 


n 

V 


n 


5 L ^ 


)) 


( 16 ) 


is average cosine 


ijli 

of scattering angle of i element 


(u-Uq) is lethargy change during slowing down = In . Val’-ie 


of ca; 
0 


•n be fixed at 2 MeV bu L value of B depends on 
temperature of the medium (E = IT) where 1- is Boltzmann 
c one t ant. 

Prom Eq. (16) we will get expression of t ' as a function 
of U/. , T, V*/ , thermal temperature of neutron. 


2 = 2.7 Modification of criticality equation : 


Major inadequacy in the four factor formiula is omission 
of non-thermal fissions. In Oklo reactors moderation was not 



47 


very effective with the result that thermalisation of neutrons 
was not complete. An important consequence of this is that 
there was significant contribution to the multiplication from 
fission in U-235 caused by absorption of epithermal neutrons. 

Oklo reactors being a slightly enriched sj'stem (enrichment 
near 37- ) this becaiies all the more iiup or tant. Here we x-fill 
develop dependence of neutron multiplication factor on epithermal 
fission of U-235. 

In order to establish a physical picture of the factors 
included, neutron cycle given in Fig. (12) is developed. Let 
at point A one thermal neutron is absorbed. As a result pf 
fission neutrons will be produced as shcvm by line AB. lo these 
fast neutrons x fission neutrons will be added by the epithermal 
fission of U-235, as indicated by line BE, so that the total 
number of fast neutrons at B is (pf + x) , As a consequence of 
fast fission effect in U-238 total number of fission neutrons 
available for slowing down is (pf + x) e. Process of slowing 
domi is arbitrarily divided into two parts, first from C to D, 
taken to be 0,625 eV , where radiation capture in U-235 is 
significant, and then from D to B, where the energy is thermal. 
The number of neutrons reaching D is equal to (pf + x) e P 20 ? 
where P^^ is the non-leakage probahility of neutrons with ener- 
gies above D and is the neutron resonance escape probability 
for uranium -238 in the smme energy range. In slowing down from 
D to E, the following additional factors are ixitroduced; Pg, 
the non-leakage probability for neutrons in the energy range D 






Fig, 12: Modi fied' neutron cycle. 




49 


D to E; P 25 > capture escape probability for V-235y and uhe 

probability that neutron while slowing down will escape captui'-e 
in various materials other than fuel. Pinally, introduction of 
P^, the non-leakage probability of thermal neutrons completes 
the cycle. The effective multiplication factor is given by 

^eff ~ ^1 ^ 2 ^3 ^28 ^28 '^’25 ^c 

Since p 2 ^ fa the capture escape p' obability for U“235i; 
the fraction of neutrons reaching the point D which undergo 
epithermal absorption in U -235 is (1 - p 2 ^) • Consequently, 
the number x of fission neutrons produced by this process is 

X = ^2^(1 “ P25^ [ (t) s P28 ^ 

where p 25 number of neutrons praluced per epithermal 

neutron absorbed in U-235. From above relation it follows 
that 

(pf+z) = pf[l - e Pp P 28 " P23^^'"^ 

Substituting it into Bq. (17) yields 

Keff = P2_ Pp ^“3 U ^ ^ P28 Pc ^ P28 

X'j’here 

-1 

P — P25^^'"^^1^ 28 ^25 (i9) 
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With P2_ ^2 ~ ®^P and P23 P23 = Pps (1®) reduces 

to 


K 


p f e p 


28 


exp (-t'B^) 


eff 


1 + 


( 20 ) 


Equation (20) is quite similar to original criticality 
equation (5) except for factor p, which tolies ca-re of epithermal 
fission and capture by U-235- Foreevaluation of p, and p^g 
have to be evaluated, p 25 is constant. p2^ and e have been 
formulated in Section 2.1.4 and 2.1.5. P-, is non-leakage 

probability for on orgy range 2 MeV to 0.625 eV and is given by 


= exp [-a-'(2MeV - 0.625 eV) ] 

where expression for t^'(Eq - E) is given in Section 2.1.7. 

Evaluation of P23J neutron resonancL. escape probability 
in energy range B to D for U-238 is diftievd-t to calculate. 

XI 

It is possible to write ppg = P28’ where n is imperically 
derived constant. It is found that best fit to the data is 
obtained by taking n = 0.5 [37]. This would imx^ly? 


. _ 0.5 

P28 ~ ^28 


( 21 ) 


where pevaluation of p2Q is elaborated in Section 2.1.4. 

How we are in a position to express s-s a function 

of uranium constaiit U/. , geological time x, poison constant x, 

volume content of water in rock pores Y'/ , content of water 

2 

structurally bound with rock, geometric buckling B , an^ 



51 


taiaperature of thennal neutron, T'anperature dependence 
of is more complicated and is dealt in next section. 


?o2 Temperature dependence of 


eff 


Major dependence of on temperature comes from 

dependence of microscopic absorption and fission cross-sections 
of elements in reaction zone on temperature as well as tempera- 
ture dependence of resonance integrals, for Oklo, expansion 
of unbound water in pores, m in equation (4), also is significant. 
Lind of formalism we have adopted in Section 1.6 does enable 
us to take into account expansion of water in pores, which 
u'-^uia be in sunercritical state. To take this into consideratio:.-. 
would require knowledge of volume of pores in Oklo rocks, V”/ 
which is not available. As we are keeping V/. as an open variahe 
this drawback d oc not inhibit us from studying criticality 
condition accurately. Study of reactor behavior after reaction 
has began will not be complete unless expansion of water is 
taken into account with rise of temperature as fission proceeds. 


2o2cl Cross-sect ion dependancie on temperature 


Absorption C3':’Oss-section of all the elements present in 
reactor zone except U-235 have been assumed to behave according 
to 1/v l"w in the thermal region. Averaged thermal cross-section 
for these elements is given by 


a (th) 


1/2 

( 295 . 6 / 1 ^) 0 ( 2 , 



library 

K^nO'ir. 
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o{2,200) ia cross-sec cion for neutron with, speed 
2200 m/sec and is given in appendix (I) for all elements of 
concern. is neutron temperature on the Kelvin scale. 

Vrc/O is a correction foi’ Maxwell distribution of thermal neutrons 


■-i:c ec 


cross- 


1 or U--255 Wsstcott correction have been used to evaluat 
ive absorptio.n and fission cross-sections. Effective 
section according to hestcott’s corn/cntion is given by 


e 


a 


eff 



+ r S(T^^)] o (2,200) 


where is a measure of the departure of cross-section from 

dopcndonce overage over Maxwell i.an canponent and can be 
■ .0 o>: j-ffa t ed ['j 8 ] as 


5^5 ^ ... 0.7245 x lO"'"^ 100) + 0.952 


= - 0.10335 X 10“^ (T^ - ICO) + 0.9495 


for absorption and fission, respectively. SCl^) is the contri- 
bution to tho total reaction rate from tho epithermal part of 
the spectrum minus Mcx^.^eliian distribution and can be approxima- 

cl c'lS 


= 0.105 X 10"'^ - 100) + Q.1185 

= 0.32 X lO""’- - 100) + 0.014 


''or absorption and fission, respectively. 
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Parameter r is used to normalise the relative contribu- 
tion of tile epithorraal and Maxwellian cODroonent of neutron 
spectrum. Value oi r was estimated by irradiating different 
samples of Oklo rock in Triton reactor by French, study group [3 S 
to be bet-ween 0.10 to 0.25. 

^n neutron temperature. In Oklo neutron temperature 
is more than 150°0 in excess of therma,! temperature of melieu. 
This is because of iiigh thermsil a.bsorptiOii ;in,d incomplete 
moderation. Relation botvreen temperature of melieu and neutron 
temperature is 

1 mass number < 25 
25 K masfj number 


S, (K T ) is total macroscopic absorption cross-section 
a, m n^ ■ 

at melieu temperature = 7 ^ cr (E. T^) N. 

iTl 

T is total epithermal scattering cross-section 

n 

' ■r ; — ■‘"• T 

= Z- .K. . 

4—1 1 s/e 

1=1 

and I is average logarithmic decriment define in 


= Tj,, (1 + 0.46 ) 

(1 + 0.30 ) 

whore 

4 (K T ) 

C.t Ui 


Section 2,1.4. 


I » 
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2o2o2 Resonance integral temperature dependence 

In the formulation develou ed above neutron interaction 
process has been idealised as neutron impinging upon a scatte- 
ring or absorbing atom which is at rest. Target atom in 
actuality has a thermal motion. With respect to the average 
position of the target atom neutron will seem to have spread 
out energy distribution. It can be sho>m for a Maitwellian 
distribution of target nucleus velocities ana 1/v cross-section, 
that Doppler effect makes no difference in the effective cross- 
section of the interaction. Only for resonances this effect 
becciBes important. It lowers Breit-Wigner resonance peak 
heigl:it and spreads it but area under peak ic not changed, so 
that for dilute elements or elements having small El there 
is negligible change in resonance integral v/ith temperature. 

However for U-238j: widening of the region of strong capture 
would appreciably iC'Crcase absorption in skirts of resonatices. 

Prom experimentally arrived values of resonance integi’s! for 
U-238 with temperature, given in MMj 5800 [38] following 
analytic form is constructed 

leff(T) = 7.46166x10“^ ^ exp[0.46 log(-Y§^)+l. 123] 

where is in °C and is temperature of milieu- 

How formulation of for study of natural reactors, 

as function of uraniumj content U/. , loose water conteat V/. , 

2 

poison content x, geological time m, gecmetric buckling B , and 
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temperature of milieii is complete. In Appendix (II) is given 
constructed computer programme for study of critica-lity 
conditions. 

Variation of above given six parameters, within cons- 
traints set in chapter one, can be investigated to find out 
situation when is one or more. Effect of temperature on 

four factors and can be studied. Also other features of 

natural reactors for comparative study of reactor systems 
can be evaluated from constructed model. 

AoJ Results : 

■-.'njol Paramcitric. de:jcndGnc (3 of Neutron te'r.poratiirG “ 

In Fig. (l6) is plotted variation of neutron temperature 
with temperature of milieu for different uranium content. For 
milieu temperature of 200°C neutron tanperature is 270°C for 
U/. = 24/( . The difference is large enougli indicating that 
distinction between milieu temperature and neutron temperature 
is crusial for modelling at Oklo reactors. This difference is 
high because of large thermal absorption cross-section and less 
slov^ing down power. 

With increase in uranium content neutron temperature 
rises because milieu content decreases as is expected from 
Eq. (2). This leads to less water content as is evident from 
E'q. (4). Net effect is worsening of slowing down power of 
materials in Oklo rock leading to hi#i degree of non-moderation. 




Temperature of Melieu.Tm in °C 


Fig. 16 Variation of Neutron temperature witi 
Melieu temperature 
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Another significant parameter which dictates neutron temperature 

is water content itr.olf. For lov^er milieu temperature the 

effect of change in V /• is r.ot much hut for higher temperatures 

neutron temperature becomes very sensitively dependent on V/, 

1 f' 

It also has a dependence on B equivaleJit poison contentj n 
as strong thermal capture implies large gap between milieu 
temperature and neutron temperature heCcPi.T.e of quicker disappea- 
rance of til ermalised neutrons. Understanding of neutron tempera -• 
ture behavior is iiuportant for understanding pareimctric dependence 
of factors idiich cparacterise behavior of neutrons in reactor - 

hi? Para.i‘i-'t r Lc 'hppendena^ of ferj. ironic, fac core 

Table {'3) gi'^'OS variation of different factors which 
characterise behavJ or of neutrons \irith diiferent parameters. 

With rise in water con tent, y, neutrons are slowed down more 
efficiently. Age, t: ' , and contribution of epithermal fissicn 
of U-235y p, decreases but resonance escape probability inci eases 
Also neutrons are brought dom below fact fission threshold 
nuickly, decroasing fast fission factor, e. In thermal region 
parasitic canturo inci'eases owing to wh' ch thermal utilisation, 
f, and diffusion area, L^, decreases. For more content of water 
neutron temperature: is less and temperature dependence of fission 
and absoi'ption cross-sections of two isotopes of uranium is such 
that p decroa-ses slightly. For low temperature fission cross- 
section of U-235 decreases faster with temperature than aggregate 
of macroscopic absorption cross-section of fuel. Decrease in 
T' and increase in p and p increases but decrease in f, 
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and e decrease it. Owing to these push and pull it is maxiinum 
at seme optimal point and decrease on both sides. Content of 
water, y, is so less that this maxima is not reached in region 
of interest hence shows an monotonic increase with water 


content rise. 

With increase in uranium content, V'/. , content of milieu 
decreases implying loss parasitic capture in thermal regioj^ 
increasing f and 1 In epithermal region less milieu content 


m;.ike slowing dov/n inefficient thus incroasing t'. Rise in U" 238 
content decreases p but increases s slightly. More 11-235 conxent 
increases p. p increases with increase in uranium conient 
because of rise of .lou'! ron temperature, increases with 

rise in U/. but for considerably hi^i b/. it decreases owing 


to loss slowing down. because of decrease in milieu and water 
content which eventually makes decrease in p significant. 


Effect of change in gcolog;ical time is not 


shown in Tabic 


(3). 


With increase in geological time cnriclmnent is more 
increasing p. Rise in U-235 content would lead to more thermal 
neutron absorption decreasing but f would increase as relative 
increase in fuel absorption would more than increase in total 
absorption in themial region. Increase in geological time 
decrease aiilieu coiitent decreasing water content thus lowering 
resonance escape probability as neutron flux in epithennbJ- region 
will increase. Because of non-moderation e also increases as pull- 
ing down of neutrons below fast fission threshold becomes 
inefficient. Non-noci. oration would increase t' as well. Net 
elfoct of is L.ocause of increase in fuel enrichment which 


increases it. 
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Etfect of chan; e in poison content is also not shown 
in Table (3). It will only effect therraal factors f and 
Lowering them both with increase in poison content, 

2 

B is a geaiietric property and not intrinsic nuclear 

pronerty. It has no eifect on different factors. Increase 
2 

in B would increase leakage of neutrons out of reaction site. 

Thus will lower K 

eii 

U-'pondance of noutronio factors on temperature of tnili iu • 

Table (4) illustrates temperature dependence of different 

factors of interest for U*/ = 24*/ , V*/ = 10/ , m' = 0.09, 

—6 

X = 50 X 10~ , T = 1.8 by and d^ = 3. Relative content of water, 

y, decreases because cf expansion of water in pores. Increase 

in temperature decreases absorption cross-section of materials in 

reaction site. This leads to decrease in parasitic capture in 

2 

therraal range increasing f and I . Decrea; e in slowing down 
power because of decrease in y explains ii'Orease in e, % and p 
and decrease in p. Behavior of p is puzzling, it increases 
I'caching maximum around 340 °C and then deci-eases. This is 
because of non 1/v temperature dependence of fission and absorp- 
tion cross-sections of both isotopes of uranium. But for our 
range of interest from 150°C to 300 °C it increases, ^gff 3-lso 
increases till 240 °C and then decreases. This is a significant 
result contrary to moat of ther.nal man-made reactors where 
temperature coeffici.nt of reactivity is invariably kept negative. 
Distinct feature of Oiao reactors is that there is lot of parasitic 
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capture in thermal range which has 1/VT dependence over 
temperature. So at lower temperatures f quickly increases 
and corresponding increase in which negatively contributes 
^eff^ ineffective because of large geometric buckling 

p 

and small absolute value of 1 . Because of small change in 
quantity of water, y <j Si t lower temperature decrease in resonance 


escape probability is not much to counter increase in f. As 
the temperature is incor'eased decrease in p eventually matches 


increase in f thus br 
points to note at thi 


inging down eventually. Another 

s stage is that temperature dependence of 


resonance integral of materials other than cranium was not 


consider owing to siaall values of there resonance integrals. 

In actuality there resonance integral will increase with 
temperature, increasing rate of decrease of p but the relative 
significance of contilbution of rest of materials to p is much 
less than that of U“-238 whose temperature dependence was 
taken into account. 


2u/'r,4 Neutron balance table : 


lo aid establish physical picture of chain reaction for 

natural reactors of Oklo, neutron balance sheet is made and is 

given in Table (5). Parameters were chosen to give = 1. 

2 “P 

B was taken 10 m " for reactors of siae of reactor zone 1 and 2. 
Uy. was taken 24'/ which is about average content of reactor 
zone 1. Content of chemically bound water was taken to be 9“/, 
content of loose water 10/ , poison content 50 ppm and tempera- 
ture wh@i reaction began to be 200 °C. Time when is one was 
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TABLE 5 


NEUTRON BALMCE SHE.ET 


! THEMAL ! 

i NEUTRONS fr 

; LEAK OUT OE ' 

REACTOR . 


; 10000 
;THEEr4AL NEUTRONS i 
ABSORBED T 
' IN REACTOR ; 


I 1917 
; PARASITIC 
4 CAPTURE OF 
! THEHIAL NEUTRONS 


‘ 10009 . 

I NEUTRONS 

I THERMALISE IN’ 
REACTORS 


8083 


MELIEU 


THEREAAI NEUTRONS 
ABSOJUSED IN 
FUEL 


: 96 ' 
.NEUTRONS leak 
I OUT BETWEM 
: RESONANCE AND 
i THERMAL ENERGIES 


/■■■■' 419 

■/neutrons ADDED' 
BY EPITHEMAL/ 
FISSION OF / 
V u-235 y 


2093 

NON-PRODUO TTVE 
CAPTURE 01 
THEEtlAL NEUTf^'-LS 


±ji 

FUEL 


9696 

NEUTRON ESCAPE i 
RESONAN CE C AP TURE , 


6030 i 

••THERMAL NEUTRONS i 


LEADING TO 
FISSION 


! J«14 S 

' NEUTRONS i 

, CAPTURED IN 
■ RESONANCE ! 

i lusGioN : 


15^0 

< NEUTRONS ■E!UTER 
i RESONiuWE 
i REGION 


...... 

FAST N.iUTRONS 
PRODUCED BY 
FISSION 


: i'],[';UiRU.iS LEAK! 

1 OUT BEFORE !. 

! i'ErT'ER,! NG i 

! RESONAI'JCE i 

! REGION ■ I 


15334 I 

FAST NEUTRONS ! 
BY thermal and*’' 
FAST FISSION ^ 


--- / ■ 440 

ADDED BY 
FAST-FISSION. 
OF 

U-238 / 
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found out to be 1741.22 million years and path of neutron 
cycle was followed to give this balance table. With the help 
of this table one can comparitively study position of nature.! 
reactors with respect to man-made reactors, 

Heutron prod-uced by fast fission of Tj-238 is 5/. of 
neistron produced by thermal fission of 0-235. 4.73*/. of fast 

neutrons are produced by U-235 epithermal fission. 15/ of 
thermal neutrons are captured by milieu and 20.5/. are non- 
productively captured by fuel to produce U-236 and Pu-239. 

Only 0.8/. neutrons 1 e&k out of reactor f i om the time fast 
neutron is born to the time it is absorbed in thermal region 
after slowing dovra. Parasitic capture ol thermal as well as 
enithermal neutrons is significantly higji in Oklo reactors as 
compared to reactors made under intended care of man. Leakage 
of neutrons is extremely small compared to modern sli^tly 
enriched reactors where 5 to 15/ of neuti’ons leak out. 

Table (6) lists some reactor physics features of Oklo 
reactors foi* same ccinbination of parameters as in neutron balance 
sheet* Average logarithmic energy decrement is quite near tnat of 
I).J) which is 0.57. 'Also number of collisions required to slow 
down neutrons is close to ^ 2 ^ which is 3b* Slowing do^,^m. time 
for D ^0 moderator is 0.5 x 10 ^ sec for Oklo it is 0.476x10 sec. 
But moderating power is much hi^er than that of pure which 
is 0.18, This is bccaijise total number density of elements in 
Oklo reactors is 0 . 2 x lO^'^ where as for ' 1 ) 2 ^ of density 1.1 gn/cm 
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TABLE 6 


SOME RE AC TOE PHI SICS EEATURES OF OICLO REACTORS 


Average logaritlimic oaergy decrement, | = 

i'lumber of collisions before neutron slows down = 

M od era tin iC P ow er , Iff 

Moderating ratio, 52;^ / - 

Enrichment of iuel, e ' = 

Total number densi'i;.j' of uranium, Nu = 

Total nmber density of all materials = 

th 

Macroscopic thermal absorption cross-section, = 

Macroscopic epithermal scattering cross-section, z-g = 
Thermal Macrosc opic transport cross-section = 

Thermal diffusi onilength, V'L^ = 

Slowing down length, Vt' = 

Migration leng!>th, V(Ii^+t') = 

Thermal prompt neution life time = 

Slowing down time = 


0.5643 

32 

0.59 

3.79 

2.9283*/ 

0.625x10^^ 

0. 2x10^^ 

0.15567 cm‘”^ 

1.04562 cm"'^ 

70.993 cm"^ 

0.9781 cm 

3.4 cm 

3. 535 cm 

0.681 4x10“ ^s 
sec. 

0.47565x10“'^ 
Sec . 


m' = 0.09, = 3 £aa/cc, V/ = 10/. , U/. = 24'/ , T^ ='200°C, 

= 10 m“^, X = 50 ppm, r = 1.74122 billion years. 
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it is 0.0334 X 10^'^. Because of high number density in Okie 
rock slowing down length is only 3.535 cm where as even for 
pure water moderator it is 5.568 cm. Therms,! diffusion length 
is considerably small because of hi,^ thermal macroscopic trans- 
port cross-section. Thermal macroscopic absorption cross-section 
for Okie is very close to that of Iodine of density 4.93 gm/cm 
which has macroscopic absorption cross-section O.I 64 cm”^. 

Thermal prompt neutron life time is less than that of heavy water 
but more than that of Beryllium. 

2o3.5 Study of criticality conditions : 

Now we go over to reflect on parametric conditions undei- 
which Oklo reactors could have become critical. 

In Pig. (17a, b, c,,d) is plotted as a function ot 

geological time, t, for different volune fraction of loose ’■Mter 
in rock pores, a, b, c, d graphs are for poison content of 

0 ppm, 50 ppm, 100 ppm and 150 ppm, respectively. Sadimentary 
assonblage of 24% uraaiima by weight of the size of reaction zone 

1 and 2 could have becesne critical as late as 1.67 byea,rs 

ago there i^as no poison and loose water occupied 12*/. volume 
of rock- If poison content was equivalent to 15 O ppm of B^^ and 
V /. = 6% then reactors could 'nt have worked after 1.94 billion 
years back. At 200 °C of milieu temperature and m' = 0.09 in rock 
of density 3 gn/cm^ possible range of time in which Oklo reactors 
could have worked by latest is 1; 67 billion years. At any 
situation in this rar,i.ge assemblage could have become critical 






in Toelow given period 


for 

no 

p oison 

t 2 

for 

50 


T 0 

f or 

100 

ppm 

T 0 

for 

150 

ppm 

% t 


1.67 to 1.795 'billion years 
1.685 to 1,81 billion years 
1.7 to 1.85 billion years 
1.716 to 1.845 billion years 


Increase in aiilieu temperature woiild shift these curves 
slii^tly to the left side of graph till a'.;out 240°C and then 
again shift them to right at slightly lovjer rate if temperature 
was still higher. Change in limits of periods would for lowest 
temperature 160 '^C shift about 0.015 billion years and fox* 
optimal temperature of about 240 °C it will shift left by about 
“0.012 billion years and then again boundaries would shift 
towards right aquiring almost above plotted position if tempera- 
ture of milieu is around 300°C. t is not very sensitive to 
temp erature. 

N'otisable feature of Fig. (17) is sensitivity of geolo- 
gical time to water content. In Fig. (18) is plotted with 

net water content per imit UO 2 S for dili'erent uranium content 
U*/. , at 1.8 billion years ago. In making this plot Y was taken 
as an open variable disregarding fonualism for it developed 
in Section 1. 6. 

For y = 0.5 assemblage could have beceme critical for 
even 16'/. uranium concentration by weight at time 1.8 billion 
years ago. It would have been critical with lesser uranium 
concentration before 1.8 billion years ap.o. All the six fossiles 
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enearthed at Oklo have average relative weight content of 
nranium more than 20/i . If reaction began 1.8 billion years 
ago with U*/ = 16°/ and y = 0.5 then how they have such high 
content can only be explained by assuming 4/ to 12/ more 
deposition took place after reaction begaj'i so that at the 
closing of chain reaction reactors had 20/ to 28/ uranium. 

If no significant deposition of uranium occur ed then for 
T = 1.8 billion years, y should be less than 0.3- for v 
greater than 1.8 billion years y would be even lessei''. Hence 
to satisfy criticality requirement of Oklo V/need not be 
greater than 12/ for d^ = 3 (from fig. (15c and 15d) .). So 
there was no shortage of water at Oklo when the reaction had 
to began only if our assumption that volume content of pores 
was between 8/ to 12/ is right. It is quite likely that 
Oklo rock lying at the base of main slippage dome in precam" 
brian B orogeny had developed many microfractures to hold 
8/ to 12/ water. 

for smaller vclues of y reactor could have began earlier 
than 1.8 billion years. As fig. (17d) indicates it would 
have at the earlicvct began by 1.85 billion years even if poi.son 
content was very high 150 ppm. If parameter temperature is 
increased it would pull up curves of fig. (18) in the beginning 
but as temperature is increased more these curves will be 
pushed down. As does not have monolithic dependence on 

temperature the effect of temperature would not be very pro- 


n ounc ed . 
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GeoniGtric buckling and criticality : 


Size of reaction zones 4» 5 6 are smaller than 

that of zones 1 and 2. Hence their geometric buclcling is 
larger than that ol zones 1 and 2, which is around 10 m“^. 
can not be assigned any fixed value as there is no certainty 
as to whether reaction began in whole area of reactor fossile. 

It could have very well began in smaller volume of concentrated 

uranium. This would mean geometric buckling when reaction began 

_2 

could have been larger than 10 m even in reaction zones 1 
and 2. 


In Fig. (19) is a plot of geometric buckling, with 

uranium percentage, U/. , for different geological time for a 
case when = 1. Condition = 1 when imbedded in 

criticality equation (20) transforms equation (20) into an 
transcendental equation for different eters. This 

transcendental equation was solved by siiiple itteration scheme 
so that we can study more accurately conditions under which 
reaction of Oklo began. Fig. (19) and Fig. (20) were drawn 
by such computation. 


Fig. (20) gives of with UX for different 

geological time, t, for = 1. ¥e had mentioned before 

that all the conclusions are very sensitive to loose water 
content in pores of rock which is again illustrated in Fig. (20). 
For very less T). reaction in 24/. uranium assemblage would be 
possible near 1.9 billion years. As summarised in Table (1) 



uudd os = 
OoOOZ = 






Fig. 20 Variation of Geomatric Buckling with 
loose water content 



75 


01 ae,o of "rimary mineralisation was suggested to be 
bo. 11 1.9 billion years. Then it must ha.ve taken couple of 
teii3 01 million years for sediments to reach such depth that 
d Gsilicanisat i on va?.s nossible which was esfjential for reconcen- 
tration of uranium -Lo take certain privileged aone to critica™ 
lity. Also conclurion from Pig. (iTa^hjC’sd) suggest age of 
start oi reaction could not have been before 1.84 Mllion ;years. 
In this li^,^t it looks likely from Fig. (20) that porosity 
fiaction of rock must be 8/. or more as noison content could ’nt 
have been very much less than 50 ppm. This sets lower limit 
to porosity. Upper ].imit to loose water content was infered 
from Fig. (15c and d) and Fig. (18) to be 12/ . So our assump- 
tion that loose .water content was between 8X to 127. f® 
substantiated as necessary condition for p.ossibility of reaction. 
Volume fraction lower than 8/. is not compatible for establish- 
ment of criticality and greater than 12/. is unnecessary. 

I'igure (19) gives plot of with respect to IT/, for . 
different geological time for V = 107 and poison 50 ppm. This 
graph can be seen iji the light of desilicanisation process in 
which uranium begtin concentrating. Desilicaniso,tion process 
con be seen as motion along x— axis. Fi'cm present uranium 
content of reactor fossiles one can infer than non of reactor 
began operating before reaching uranium concentration of about 
207 by weight during desilicanisation. Limit 20% is set 03/ 
accepting that upto 47 uranium could have got deposited after 
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reaction began, ^his would imply from Pig: (19) that biggest 
reaction could 'nt i'lave began some 1.79 billion years ago for 
"V ~ 10/. . Por hi^or V this upper limit stnfts down and for 
lower V this limit shifts up. Por reaction zone 3, 4, 5 and 6 
B was considerably higher upto 20 m ^ and these reactors could 
have worked around 1.79 billion years. Only when U/ reached 
around 24/ » P or a situation when all the six reactors worked 
simultcin eously we can conclude that either U/ in smaller 
reaction zones should have been more than that of largest by 

4/ or that loose water content was more in smaller zones than 
in bigger, 

Por V = 10/ smaller reactor zones with U = 24/ could 
have only become ci-itical before 1.78 billion years ago. This 
sets lower limit to time of their working. But it V/= 12/ 

(from Pig. (20)) this lower limit shifts to about 1,76 billion 
years. If reactors were working simultai'i eously than the period 
in which they worked is between 1.81 to 1.74 billion years 
with lesser water content in larger zones and more water content 
in smaller zone. 

If process of desilicanisaticn occured in very short 
time span compared to geological time considered than bigger 
reaction zones could be critical as late as 1.74 billion years 
for U/ = 23.6/ (from Pig. (19)). And for V - 12/ , we can say 
from Pig, (20) that bigger zones could have heccme critical as 
late as 1.68 billioh years. This gives another picture of 
smaller zones becoming critical first aror.nd 1.81 billion years 



to 1.76 billion years and later larger acres becoming critical 

•upto as late as 1.68 billion years with Y = 12*/. and poison 

2 

50 ppm. If we i;ake 13 for intermediate zones to be 15 
then range for such reactors become 1.81 billion j ears to 
1.73 billion years, v/e can summarizes 


f or 



10 

m“^ 

T 

1.81 

to 

1. 68 

billion 

y ears 

15 

m-2 


1.81 

to 

1.73 

billion 

y ears 

20 

m-2 

X 1 

1. 81 

to 

1.76 

billion 

y ears 


for loose water contait V from 8'/. to 12/ 
and for uraniim content 20/ to 28/ 
and tor poison content of about 50 ppm. 


For higher poison contents these limits shrink down. 
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CHAPTER - III 

STUDY OF SUSTENANCE OF REACTION 

3^1 Dynamic model of reactors 

In tills chapter a mathematical model of generation 
and burnup of different isotopes, when reaction at 

T T -'-vi o n osir Q 

has been made. This study has been motivoued by 

, _ . _ ppqctors. 

investigate sustor'ance of chain reaction xn UKio n-c 

1 -P"i 1 'i cT r & as ^ d 1 

In the below developed formalism we wonld seep imei 

as flux as open variables along with eiohi open -parameters 

« H d 1 ti oiibI 

considered in criticality model. Fo.llouu'rig set oi 
assmptions have been made 

^ com pan* ed 

1. Time for which reactors operatoo .is negl^S^ 

, . ^ -jf < 5 BC t or 

to geological time. Hence geological time a 

time have been treated as unlinked -variables. ^ 

It Or 

2, Heactofs op orated continuosuly irit'hout brea 

T j of majol 

breaks were small enough so, that naiurai aeoos) 
isotopes Coin be neglected. 

;•). Flux remain, d constant all throiigout reactor opera 


= lol Burnup/hv i Idup aqua t ions : 

^ TT and jnitiadly 

Once chair;, loaction begins burning up of u atx 

ent poisons effect process various urar ' 

plutonium isoto-es are formed, also formc-d are poisonous 

c i 3 /b^ 1 ^ 

fission products. Isotopes which are fc-:.aed in appre 
quantities such that they effect are 13-236, Poi 

Pu-*24C, Pu-241, Xe~135j Sm-'149 and other fission prod 
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Out of them PU"2:'3S' aud Pu— 241 are fi s&ioiiatbl e and. contribute 
positively to 

Distinct feature of following isotope density equatio 
is thatfluence (0) as well as flux (cp) , occurs in them in 
contrast to man made reactor burnup /buildup equation where 
only fluei-ice is a variable. I'his is done to include the possi- 
bility that reaction operation period could be comparable to 
isotope decay half l:lves. Such a possibility has been suggestei 
by the study of isotopic abundance in Ohio samples. In fact 
the lower limit to reaction operation c.eriod had been suggeste< 
to be 0.2 million years by the same study. 

Burnup of li~235 can be describod by following equation 

25 25 

d~ = “■ '^‘'25 ""a + ^29 ^ 29 /'^' " ^ 

( ^5 1^25 ^9 ^'^29 ^ ^22) 

0 is flucicco, 9 is flux, is decay constant of 

Pu-239, and sire neutron yield per fission of U-235 nnd 

Pu- 239 . First term takes care of loss of U-235 by thermal 
neutron absorption , second term gives gain in U-235 population 
because of Pu-239 decay. Third term ropresents loss of U-235 

f 

in epithermal regl-^Jir. Pj_ is non-lefikrige probability till 
180 Kev energy. Contribution of U-235 loss from neutrons born 
fron Fu-239 fiscio! is also considered in last term. Matural 



decay of U-235 is ignored as its half life is of order of 
riillion of years where as reaction period order is most 
likely to be inillion hears, following ocoetion takes care 
of Pu“239 generation in Reactors. 


dN 


d9 


~ = ^'^'28 ’ ^ V'25^f^ ■*' 


^29*^6/ “ *^29^' 29 A’ 


(23) 


The first two terms accounts for thermo,! and er r thoimal 
absorption in U--238 which leads to formation of Pu-239. Third 
term is loss of Pi'i“239 by absorption of thermal neutrons 
leading to fission or formation of Pu-240, The lo-st term is 
isotope decay term \.rhich leads to loss of Pu-239 but gain in 
U-235 density. ;[lpithcrmal and fast absorption by Pu-239 has 
boon neglected. 


Equrdion for Pu-240 is 


dR 


20 


dG 


N 


.29 


29 




.20 


20 


^20 ^^20 

Q) 


(24) 


First term in above equation is ;m capture component 
of tiiird term in oquation (23) o.nd represents formation of 
Pu-240 by therm.il rioutron absorption by Pu-239. Second torm 
gives loss because of thermal neutron absorption by Pu-240 
which leads to formation of Pu-241 as fission cross section 



of Pu-240 is negl-!-t,il)io. Last term is loss Lccause of 
Isotopic decay of F 0-240. Here also cpitLermal/fast absorp- 
tion is neglected. • 


Equation for Pu-241 is 


dH 


21 


d9 


af - H,, aP 

cl 21 ci 


‘21 ^’21 


(25) 


Por every theiiaal neutron absorbed iii Pu-240 one Pu-241 
atan is born is ^;b■.~t tlio first term Second term gives 

loss because of tboim'al absorption of neutron, part of these 
thermal absorption leads to fission. Tbird term is nat-ural 
decay term. Hero oigain epithermal/fast effect havo been 
neglected. 

U-236 is olso formed in largo quantities and its numbei 
density is constr. -ired by following oquatrcn. 


div 


_2.6 

dO 


,25 






■'‘26 


o 


26 


( 26 ) 


a 


First torn is formation of U-236 from thermal neutron 
capture in U-235. Second term is burningup of IJ--236. U-236 

fom'iatioi.1 by epiti'crnal capture in U-235 bas been neglected. 

Most important burnable poisons are small quantititjo 
of rare earths cuid other elements with strong capture crosa- 
socticn. Their presence had been idealised by B equivalent 


me- sure in Section 
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Equation for number density is 

B 


(27) 


Only term present represent burn-up of poisons. 


Eission of various isotopes (U-235!, Pu~239, Pu-241) 
leads to formation of fission products which, adds to parasitic 
capture of neutrons. Most impoi'tant isotope build-up frcxn 
capture point of view is of Xe-135 and >jm-149. Xo-135 is born 

in fission process as well as is formed as decay product of 
other fission products. Below given diagrau illustrates isotop( 
decay chain leading to formation of Xe-135, 


To 


135 

-■71 




fission 


<1 min. -pl35 “S.? hrs. 9„ 2 hrs. 

~ Xe 


fission 


;ission 


n 


155 2x10 yrs.,^, g^,135 ( 


stabl e) 


¥e can assume that all Te-135 formed immediately decays 
into Ir-135 and does not absorb neutrons to create 136 mass 
particle. To evaluate Xe-135 number density 1-135 number 
density has to be knoxm. Below given equation describes number 
density of I-135 isotope. 

dN 

= -1*1 “a I*i /-P+(Y5 °f'*25+V9°fs29) 

( 28 ) 



First and second term describe depletion of ^ 

absorption and decay respectively. In last term ar.d Yo 
are fission yields of Iodine and Tellirium combinded, from 
U-k35 and Pn-239 fission respectively. They have been multiplied 
by macroscopic fission cross-section of U-235 and Plutonium- 239. 
Once NI is kno-s./n from above equation it can be imbedded in 
Xe -135 equation given below to get number density of Xe-135. 


dJi-. 
z. e 

d0 


Ji. B cl 




Xe Ke' 


X I 






(29) 


First and second terms are decay and capture term of 
Xe- 135 . Third term is production of Xe-135 by decay of I- 135 , 
last term is fission yield of U— 135 and Pu— 239? respectively. 

Sm -149 is stable fission product. It is not only formed 
as fission product a,lone but also is a product decay of Pffi -149 
is 


.Pm 

”dG 




25 


(50) 

Pm -149 is not born in Pu-239 fission, hence third term 
only contains yield by U-235 fission. First and second terms 
are decay and neutron capture terms. 
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Bq-uation for Sifl-149 is 


d]\T 

Sm _ S) 
de “ " a 


N., 

8m 


+ 


> - /cp + Y 

Pm '' 


Sm 




(31) 


Sm'^149 is not born in U-235 fission;, hence last terra 
only contains contr i.bntion from p-u-239 fiosion. First term 
represents neutron capture by Sm-149 leadin':, to its depletion. 
Second term describes contribution from p decay of ?m~149. 

Pest of the fission products can be considered to be 
equivalent to a hypothetical element with microscopic cross- 
section of 80 barns till 10^^ fluence and 66 barns after that, 
fhe equation governing number density of fission products can 
be written as 


Cii.''! jy .-P* 

d© a 


.!■ • 'J j 


+ 2(y^*'-P*c72%2^ + Yq*^*c^,?%oq + 


29 


21 ' 


(32) 


In equations (28, 29, 30, 3l) contribution to the number 
density from fission of Pu~241 had been neglected. In equation 
(32) we have parti^.'.lly made up for it by including contribution 
from Pu-241 fission and taking cummulation yield of it’s fission 
products to be one. The first term in above equation is an 
absorption term of an hypothetically unified fission products. 
Contribution from fo.st fission of U-238 f.md epithermal fission 
of other fissile isotopes have been ne^^ected. 



I'J umber density of U~238 during reactor working can 
be described by following equation. 

pg pQ _ oc 

d©” ~ ~‘^c'^'28 ” ^ ^25 

+ ( 33 ) 

Here first tei^m represents loss of U-238 by capture of 
thennal neutron. Pirst term in square brs.ct.et gives loss 
because of epithermn.l capture and second term gives loss because 
of fast fission of U-238- 

There are six other elements of importance present in ; 

reactors. They arc Si02, AlpO^y PeO, MgOg, MnO. They will 

get depleted by capture of thermal as well as epithermal neutrons,, 
thou^ capture of epithermal neutron is quite less caapared to 
thermal neutrons as is explicit from last chapter. He^ecting 
epithermal capture tneir number densities CrUi be written as 

oxp(- cT © ) 

where N° is number density of each of them at the start of 

n ~PA PI 2 P 

reaction. As o is of order of 10 ^ fluence till 10 - 10 

a 

will make negligible effect on their number densities. Hence 
their depletion can. bo neglected if reactor behavior till about 
10^^ - 10^^ fluence have to be studied. 
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Now we have twelve first order eq_uation(22 to 55) » some of thuia 
coupled, , describing build-up and b\jrn-up of crusial isotopes 
frcm criticality point of view. 

Temperature dependence of cross-sections ‘ 

Cross sections and other data needed for solving above 
set of twelve equations is given in Appendix (I) for neutron 
'elocity of 2,200 m/sec. The temperature dependence of needed 
:ross-sections which do not obey lA law are as follows; 

u-239s 

= 1031.1(0.913+0.427 Tn/200) [293. 4/(Tn+273) ]^’^/1.128 barns 

0^9 ^ 747.73(0.924+0.322 Tn/200) [293.4/(Tn+273) ]°*^A.128 barns 

u-240: 

0^0 295(1.017+0.065 Tn/200) [293.4/(Tn+273)]®*^/l.l28 barns 

a. 

e-135: 

oj® = [3. 9753-(0. 445/400) (Tn+ 273)] 10° barns 

a 

1-149: 

= [4.2(1.304+0.289 Tn/200) [293.4/(Tn+273) ]°^^/1.128]lo\.arne 
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Tn is neutron temperature in above equations. For rest of cross- 
section ’Seeded, l/Vi dependence on temperature was used. For 
hypothetically unified fission produts, absorption cross-section 
is a function of fluence and as no analytic form for it was 
available following descrete function can be useds 


f .p, 


0 ^ 

^^/1.128 

for 

0 < G c 

°*^/l,123 

for 

10^^ C 


lo3 ‘S'c t e!7iG for numer ical soluti on : 


Proai results of Chapter 2 we imow conditions under which | 
chain reaction would have began in Oklo. Once chain reaction 
began behavior of j/eactor can be Isiown by embedding solution of 
above given set of twelve differential equations in critic.ality 
programme de’V’eloiod in chapter 2. Analytic solution of above ; 
given set of equations is pr olblimatic. Scheme for numerical ; 
solution of abc^’e equations was developed using Eunge-Kutta method 
vrhich is an fourth order approximation me'lhod. Computer i 

programme made for solution, of these equations is given in 
Appendix (III) and Wcas designed to be used c-s an subroutine for 
main programme of chapter 2. 

As the soloition of above given equations is needed over 

21 

a very large range of fluence from sero to 10 and order of 
different coefficients in these equatiens rlso is scattered very 
much lot of computational problems came up. Choise of step length 
became very important to minimise error. Dicisive advantage of 



Runge-Kutta method is that control of step length from within 
in possible so that error can be kept within desired limit. In 
Rnnge-Kutta method at every step four auxiliary quantities are 
canputed for every equation. If there are m coupled equa- 


(i) 


tions like y^'"^ = f(6, y^^^ , y^ . .y'’’'" ^ where to y 

are m var'iables, than at step of step length h, for i'^^ 
equation these four quantities are 


.O’a) 


( 1 ) 


(m) 









1 


1 .1 


“l(®n + 3 >"’ 4 + 2 \ V ^2 


11 


n ■ 


m, 1 ,m<, 

• ,yy^+ 2 


+ -2 bJ, 


'n 


+ 


1 r>i 


n ■ 


m 






= hf . (e + h, y^ + cj, . . . , y^ + c\ . . . , y^ + c“) 


"t]! 

And the new va.lue of i variable is 


^n+l 




n' 


In this schcric other than round-up error and truncation 
error, an error is caused by the fact that slope f is evaluated 


at (0. 


n’ 


1 

'n^ 


. vp instead of at (0^, y^(G^, ,y“(©j^ ) 




This last kind of error poses ccmputaticnal problcTiS, if absolute 
- lue of any partial derivative = "‘i/ p y is largo this 
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error becomes significant. In our set of equations decay 
terms lead to this problem, which is absent in set of burnup/ 
buildup equations of man-made reactors. It can be shown that 
this error could be controled if step length h is chosen such 
that quantity 


i 

n 


2 [ 


On - 


B 


1 

la 


4 


(54) 


is within, say, 0.05 ^ 0.2 for every equation. 

If any >0.2 then h must be decreased by say 50;$; and 

if ^n equaticais then h should be increased 

to say, double its value. This lower limit is to minimise 
number of steps to be computed. This ristriction on h brings 
it down increasing number of steps hence increasing roundoff 
errors. Roundoff error was minimised by working in double 

ijiii 

precision in which 16 digit is rounded of. Truncation error 
for every step in Runge-Kutta scheme is of the order of [h/G]^ 
where 6 is interval which has to be covered with steps of 
length h. Truncation error just gives upper tolerable limit, 

H, to stop length. The value of h found from criterion of 
equation (34) should in no case go above II. With all these 
criterion error per step can not go above order of [h/0]^ 

Another alternation check for accxmiulated absolute error order was 
followed. Equation for equivalent poisons has simple expo- 
nential aniolytic solution. Result' of analytic solution and 
numerical solution were matched for assessment of absolute 

rOiindnn nrid tn in no "hi on Arnnn 
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3o2 Results : 

jo2ol Fluence achieved after criticality : 

A situation was chosen when reactor becomes just critical 
and for particular assumed flux, 9, duration till hecomea 

less than one was computed using conputati onal module developed 
for burnup study. Duration of reaction, T? for different 
assumed flux, 9, is plotted in Dig. (21). Both axes are on 
logarithmic scale, fluence (= 9-1) reached by reaction was seen 
to vary from 10^"^ to 10^^ neutrons/cm^, depending on flux. 

Pluence was less for lower values of flux. Ihis is illustrated 
by the deviation of actual curve from dotted line which represents 
relation between T and 9 for constant flu once of 10 n/cm 
Reason for this deviation is that for lower flux there is more 
natural decay of isotopes, in other words there is more Xe-135 
and Sm-149 production as well as less number density of Pu-239 
and Pu-241 for same fluence. All of these changes contribute 
negatively to neutron multiplication factor so that net fluence 
reached till K beccmes less than one is less. Humber density 
of different isotopes of importance at fluence 10 n/cm for 
different values of flux is given in Table (7). Here we can see 
number density of Xo-'135, Sm-149 and Pu-241 decreases with 
increase in value of flux. Xe-135 and Sm-149 number density 
decreases faster because of smaller decay constants of I 135, 
Xe-135, Pm-149 and Sm-149. Number density of Pu isotopes show 
insignificant change in the range of flux 10 to 10 neutrons/ 
cm^-sec. In this flux range change in Pu-239 and Pu-240 was 



o 



Fig. 21 Relation of the time till which Reactor operates -.Yitho-.ii 
further accretion of Uranium with Neutron flux 
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TABLE 7 


: ELUX ; 

: 9 ' 

( n /cm^/sec ■ 


10 


10 


10 


14 


^ELUENCE -ic ^ TQ 

0 10 ; 10 ^ 

I n/cm^ ’ ! 


10 


18 


10 


19 


10 


19 





I 

4 

} 


:^''U~235 i 0.182475x10^^1 0.182475x10^^! 0.182475x10^^; 0.182475x10^^ 

! ' ‘ 


pPu- 23 S ‘ 0 . 35484 x 10 ^^ 0 . 35484 x 10 ^^ ; 0 . 35484 x 10 ^^ ! 0 . 35484 x 10 ^^ 

^ J j 


'‘^Pu--240 : 0.2608x10^^ ! 0 . 260836x10^^' 0. 260836x10^^1 0. 260836x10^^ . 

{ ^ — i '""'i'''" f ; 

! * \ 

I^Pu-241 : 0.80826x10' i 0. 80228x10^^ , 0. 8022825x10^1 0. 8022825x10^ . 


I^'^U -236 ‘ 0 . 92 x 10^5 : 0 . 92 xl 0 ^^ . 0 . 92 x 10 ^^ ; 0 .g 2 xl 0 ^^ 


’ 0.312448x10^®; 0. 312448x10^®; 0.312448x10^®; 0.312448x10-^® . 


:^^I-135 ■ 0.1987x10^^ 1 0. 246426x10^^1 0. 2479674x10^^; 0. 24796723x10^®! 


i%e-135 

,■ 0. 271x10^® 

0. 2336x10^® 

! 0.7979x10^^ 

! IP- 

; 0.79776x10^^^ ; 

Fpm-149 

t 

] 0.364x10^® 

0.56865x10^"^ 

i 0.56909x10^'^ 

• 1 

; 0.5691x10^"^ 1 

1 

:/^Sm-149 

f 

“T — — n 

! 0.532x10^^^ 



0. 453897x10^^ 

1 9 

; 0.96757x10^ 

1 ■■ ^ 

! 9 I 

; 0.96357x10^ i 

[ - , ^ . \ 

fi'.p. 

i 0.752x10^^ ^ 

^ 0.752x10^^ 

, 0.752x10^® 

1 

1 c ! 

10.752x10^° i 

1 , 1 


Uy. = 24y , X = 50 ppm, B^ = 10 = 3 gn/cm®, m’=0.09, 

vy = loy , T = 1.74122 Billion years, T^ = 200^0. 
Conversion ratio = 0,608846 is constant. 
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detected to be of oraer of 10^ and 10 respectively. This 
is expected as order of Pu decay is about (10“^^ EPu/t) , which 
is negligible compared to other terms in higher flux range but 
becomes significant at loiter flux range when reactor working 
time becomes comparable to decay constant at Plutonium. In 
Pig. (21) duration oi reaction for flux lower than 10^ n/cm^-sec 
is not plotted. Ccciputation of these points needed extreaely 
large number of itterntions to keep results in desired error limit. 
This was because coiatrolled step length of fluence was small 
owing to increased significance of decay terms. Even to evaluate 
duration of reaction for 10 n/cm flux about 0.3x10° itteratious 
were needed to keep error in desired limit. Por evaluation of 
above given results error in criticality situation chosen is 
equivalent to U/. = ( 24 + lO"^)*/. . 

j r, 2 o 2 Additional daposition needad for siistcna.nce of reaction ■ 

In last section it has been demonstrated for one likely 

situation that fluence after starting of reactor does not exceed 
17 19 2 

10 to 10 . n/cm range. Por reactors to have worked till 
fluence of 10 n/cm‘^ as suggested by isotopic studies [42]. 

Some mechanism of compensating for loss of criticality is needed. 
As vras pointed out orrlier situation of criticality was achieved 
during desilicanisation process in which uraaiium. got reconcentrated 
in reaction zones. ¥e hypothesis that after reaction began 
desilicanisation process continued and uranium kept on depositing 
in reaction zones. In this section we will give computed values 
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of quantity of additional mineralisation needed to sustain 

21 2 

reaction till fluence of 10 n/cm for different initial 
criticality situations and assumed flux. 

For criticality situation with U/. = 24^ additional 
mineralisation needed to sustain reaction has been tabulated 

p IQ 2 

in Table (8) for assumed flux of 10 n/cm"' and 10 n/cm . 

Quantity of uranium added in reactor zone is represented by 

increased U/ (content of uranium by wei-flit in percentage). 

Here we note that for higher flux less additional mineralisation 

is needed. This is because more PU“259 and Pu-241 are produced 

for higher flux leading to more fission. At fluence 10 n/cm 

20 3 7 

PU“2:59 and PU“241 number density is O.3413rl0 and 0.967x10 ' 
for flux 10^"^ n/cm^-sec and is 0.864x10^'-'' and 0.11x10^^ n/cm^ 
for flux = 10^^ m/cm^-sec. 

In Ohio reactors only in seme localised area uranium 

content was as high as 60/. [9 ]• Otherwise content of uranium 

was scattered between 20/ to 50/ . The results in Table (8) 

shows that even for Sidditional mineralisation of 6/ fluence 

achieved is near lO*"^ n/cm^. This indicates that initial 

situation of 24/ uranium could at the most be a localised 

situation in reaction zone. For four more examplary situations 

additional mineralisation needed hiive been cemputed for two 

p p 14/2 

values of the 10^ li/ciii -sec and 10 n/cm -sec and are 

in Table (9). Initiol situation with lower U/ were chosen to 
allow larger addition of uranium. Two situation with' hi^ poison 
content, x, were chcisen to allow for help in sustenance by 



TABLE 8 


U/i =24/. , V/,=10/ , X = 50 ppm, = 10 


T^ = 200^0, = 39 cm^, m^=0. 09, t’=1. 741 22 byrs. 


^ PLUX 

1 

j 

10® 

i 10^^ I 

> 

1 - - .... - . 

\ 

1 Pluence at which 

1 mineralisation 

I is resumed . 

0.7026x10^® 

' ? 

i 2.^ • 

; C. 15x10-^" 1 

1 I 

! 1 

: PIUMTCB ; 

■ - - , j 

U/. 

U/ 

{ 

j 

i 

i 

PO * 

0.1x10'^^ i 

25. 25*/ 

[, I- ^ r n i . „ i.,i i . iiiin 11,1 

) 

24.25/ 


: ' i 

0.25x10^° ! 

1 < 

26.30/ 

r " 

l 

1 

24.59/ 

j 

; ’1 

i 0.5x10^'^ i 

- - - i 

27.38/ 

' " ' "■■■'"' ” 

26. 25/ 

i 

! 

i 

t 

1 

■ 0.75x10^^ ■ 

i ...1 

' ' ‘ - "™' 

j 

28.14/ 1 

25.9 / 

1 

1 

j 

! 0.1x10^^ ! 

! i 

j - 

28.84/ i 

26,62/ 

j 

j 

■ 0. 25x10^^ 1 

; j 

"1 

1 

i 

32.90/ 1 

i 

30.82/ 

j 

i 

1 

i 

1 

’ 21 ^ 
■ 0. SxlO'"'^ 

39.04/ 1 

37.76/ 1 

1 

0.75x10^^ 1 

- - - - -J 

1 

45. 68/ 

t 

44.35/ ! 

f 

1 

! 9 9 i 

; 0.1x10^^ 

t 

f 

51.71*/. ! 

50.30/ ; 
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'burning up of initially present poisons. Also a situation was 
chosen with high poison content and smaller size of assembly 
to acccmodate possibility that reaction could have began in 
smaller area which is capable of susteinance as well as can 
feed neutrons to the rest. From Table (1) we icnow that age 
of primary mineralisation falls around 1.9 billion years. 
Considering results ol Chapter 2 we can conclude that it is 
unlikely for Oklo reactors to have become critical for U/. less 
than 16/.. Two situations for 'Q/ = 16/ and other two for 
U/ = 20/ were taken. 

In none of the situations in Tabic- (9) fluence would 
21 2 

have reached 10 n/cm within constraint on uranixma content 
infered frcm fossilo of Oklo reactors. Fluence reaching 10^° 
n/cm is possible with 2.5/ to 6/ additional mineralisation 
for initial uranium content of 15/ to 24/ , V/ = 10*/, poison 
content of 50 ppm ajad flux range 10® n/cm'^-soc to 10^^ n/cm^-sec. 

R p 

If neutron flux in reactors was less than 10 n/cra -sec than 
this rajige shifts to hi^er values and shifts to lower values 
for higher flux. Another point to note Lore is that additional 
mineralisation needed is less for highei’ Initial uranium content. 
Considerang constraint on uranium content of Oklo reactor fluence 
till 0.3x10 setans possible (in situation (c)). For higher 
poison content fluevico reachable in reactors increases. Situation 
(a) and (d) have extremely high poison content. From (a) and (d) 
we can roughly draw upper limit to fluence reached as around 
0.5x10'^'^ n/cm*^. This needs 8/ to 12/ extra mineralisation 
with initial ura.niuni content between 16/ to 20/ . We can 
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conclude that in most of the reaction area fluence between 
10^^ n/cm^ to 5x10^*^ n/cm^ would have been reached with at 
least 2.5% to 6% additional uranium for initial uranium 
content between 16% to 24% and poison content of 50 ppm or 
more. 

Seme localised areas in reaction zones have been 
reported to have had more than 30%- uranium content. In 
reaction zone 2 fairly large area had uranium concentration 
up to 40% , at some spots in it concentration is 60% . in 
these locations there is a possibility that larger value of 
fluences would have been reached. 



CHAPTER - IV 


4„I On general possibility of nuclear chain reaction in Nature : 

Apparently the only major roles of the fission process 
in nature are to pro\"ide an upper limit to mass number and to 
some extent influence the abundcnce' of fission products. Apart 
from this we may forget about fission if we are interested only 
in the major features of our world. Ihe fact that the number 
of neutrons emitted per fission, , is around two seems to be 
of very minor importance to nature. Also thei’e is no pressing 
theoretical reason that ' be greater than one. In fact one oj 
the earliest nuclear researcher Victor Wiesskopf reflecting on ; 
pre-second world wai’ days had written " Many of us hoped that 
the number of neutron per fission would be low enough to prevent 
making of a bomb" [52]. But ” is large eiiough for nuclear 
chain reaction to be powssible. How we Icnow at least in Oklo 
nature had made use of this possibility. Other possibility of 
occurance of nuclear chain reaction in nature are there. This ; 
in the period following formation of earth. One school of : 

historians of earth feel that earth was cooler when it first j 
formed because Neon is extremely rare in our atmosphere but i 
comparably abundant in stars. Neon's weight is same as moleculi 
weight of water and vrater differs from Neon chemically in that , 
it reacts under 'cool' conditions to form compounds where as i 

i 

Neon is inert. , This would mean rejection of hypothesis | 
that earth's mental is In the process of formation of I 



earth from plasma cloud most of the uraniiim xfould have been 
bracketed between iron core and surface because of gravitationa 
sorting and ionisation energy diff erenciati on. In the process 
of crystallisation of rocks under pressure uranium rich fluid 
would have been expelled [9 ]. As U— 235 enrichment at the time 
of formation of earth 4.5 billion years ago was 23'/. relativel; 
small mass of uranium would have gone critical. But fossil es 
of these reactors, if at all they exist, would be still buried 
deep under crust of earth. Unless some orogeny related phenomer 
has led to their exposure there is little possibility of detect: 
them. On surface of corth till the end of fledean period 3.4 
billion uears ago uranium remained as dissauinated UO 2 in crystr 
nine rocks as there xfl'ere no known mechainisms which could have 
concentrated uranium. During Archean period (about 3.4 to 2.6 
billion years ago) procaryotic life forms evolved but there was 
no significant amount of free oxygen during this time and uranio 
remained as un oxidised UO^. In this period uranium did conCentr 
by’' means of gravitational sorting of heavy stream sediments, whi 
included UO 2 . These unoxidised ores kno-^m as paleoplacer or 
quartz-pebble conglomerate ore rarely contained more than few 
tenths % uranium and oven thou^ relative U-235 content was 15/ 
to 8/ there was total insufficient uranium concentration for 
criticality. Then during Proterozoic (about 2.6 to 2.0 billion 
years) important changes occured. This was the period of intcjnc 
mountain building and subsequent weathering. Evolution of 
oucaryotic life capab].o of photosynthesis librated oxygen and an 
extensive surface *i‘ocervoirs ^ of un oxidised uranium bearing 
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sediments and orogeny related exposures were attacked by oxygen 
and large quantity of ccmplexed uranyl ion were hydrocharically 
released. This uranyl precipitated as uraninite, pitchblende, 
or coflinite often resulting in large massive views of uranium. 
Theyrrepresent world's richest reservoirs. Though U-235 conceE 
tration was 4/. to 2/. reconcentration of luranium because of _ 
tectonic activity could led to criticality. Oklo reactors 
evolved in this period and there is no special reason why many 
more such reactors would not have worked in that period. Next 
major r ec one en trot ion of uranium occur ed when plants evolved. 

The organic debris of plants in terrestrial sediments lea.d to 
formation of pitchblende vein deposits. But it wa,s too late. 
Concentration of U-235 -had dropped to a level when chain reactio 
was possible only in some specific geometries. Several years 
ago some petrified trees were discovered in the Colorado mine 
in the vicinity of the San Migual River deposit. One of the 
trunk v;as a^bout 30 meters long and a more than a meter thick 
and other a little less. Over a hundred tons of uranium ore v/as 
obtained from these trunks .' From this amount of ore seme six 
tons of U0o> about 2 gms of radium and number of other elements 

C. 

were extracted [53]. Had there "been trees in proterozoic period 
there would have been thousands of ilatural nuclear Reactors. 

. 2 Co nclusions : 

Effective neutron multiplication factor, ? does n 

have a monolithic dependence over temperature. Till about 
temperature of 240°C temperature coefficient of reactivity : 


( 1 ) 
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positive after that it becomes negative (from Sec. 2.3.3). 

This will have important implication for deciding flux level 
reached in reactors. Plux is dependent on excess reactivity 
wdich in case of Oklo reactors is because of rise in temperature. 
If initial temperature of milieu was more than the tarperature 
at which is maximimi then reactors would have worked at 

very low flux. 

(2) Slowing domi length and thennal diiiusion length are 
considerably s}iii.(.ll in reaction zone (from Sec. 2.3*4). If 
outside reactor zone thermal diffusion length is small, 
which is likely to be the case, \-re can say that neutrons 
fran our reactioai zone did not reach other in appreciable 
quantity. This wouldindicate that reaction zone 3-4 and 
5-6 (in Fig. 6) were not coupled reactor systems. Such 
possibility was hinted by Maurett [9]. 

(3) In Oklo j'eactors leakage was verj" less, ■ 1/. , because 

of its unique moderation feature and high parasitic capture. 
This has a massage for radiation safety concern. Average 
logarithmic energy decrement, number o± collisions needed 
for slowing doxjn and slowing down tiiue of Oklo moderator is 
sli^tly better than that of D^O moderator but moderating 
power is much bettor than even that of natural water. 

This is because of high number density of Oklo moderator. 

On the other hajid parasitic capture in thermal as well as 
epithermal region is quite high, It_ is as if Oklo moderator 
absorped neutrons which were to have leaked out. This points 



103 


to ail new reactor concept in which leakage is quite less 
at the same time it serves a purpose of production of 
useful isotopes. 

(4) following set of conclusions can be drawn about conditio 

under which Oklo reactors became critical ; 

(a) Criticality in Oklo was possible before 1.67 billion 
years. Tlu s lower limit on g'^^ological time is much 
less than age of primary mineralisation. 

(b) Volume ft'O-otion of loose water (7/. ) greater than 6/. 
is neCGSs ry for criticality. In other words volume 
fraction of microfractures in which water could stay 
should bo more than 6/. for criticality to be possible. 
Criticality condition in general v^.s found sensitive 

to water content in reactors. 

(c) Most likely reaction did'nt bv-giii before l6/. uranium 
content by wei^t had been accumulated in the process 
of d esilicnnisaution. 

(d) For all six reactors to have worked, simultaneously 
either urajiiura content (U*/. ) in smo-Uer zones should 
have been 4/. more than larger zones or loose vrater 
content sJiould have been moi'e in smaller zones. 

(e) For loose v/ater content 8/. to 12/. and uraniium 
content 20/. to 28/. and poison content of 50 ppm. 

for - 10 m*"^ x 1.63 billion years 

= 15 X > 1.73 billion years 
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p ^ P 

B = 20 m T '>1.76 11111011 years 

The upper limit Is left open because of uncertainty 
In age of primary mlner:--llsatlon. 

(5) Onco reaction begins functioning fluence between 

IV P IQ P 

lO"^ n/cm^ to lO"^" n/cm^ would have been reached without 
further accretion of uranium. 

20 

(6) In most of the reaction area fluence between 10 

P PO P 

n/cm*^ to 5x10 n/cm would have been reached with at least 

2.5X to 6/. addltlonca uranium for Initial uranium content 

between 16*/. to 24/ and poison content of 50 ppm or more. 

Fluence attained is sensitive to poisons. Maximum fluence 

reached throu^ continued mineralisa.tion considering bound 

on content of ursaiium inferred from fossiles of Ohio reactors 

is half of 10^^' n/cm^. Study of isotopic content of Th, 

21 2 

Bi, Ru and Md [Al] suggests final fluence to be 10 n/cm . 

In that case hypothesis of continued mineralisation alone 
can not explain sustenance of reaction. 

Shortcomings of this study ’ 

Major drawback of this study is that thermal and epithermal 
flux aro not differentiated. In other words epithermal flux is 
oxa.ggGrf.it ed . This lea.ds to over estimation of epithermal quantitie 
Another weak point in criticality model is calculation of fast 
fission factor for which very simplified formalism is used. 
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Another limitation of this study is that spacial distribu- 
tion of different elements in reaction zone has been ignored. 

As far as study of criticality conditions on.d reactor physics 
features of reactors are concerned this is not much of a handicap 
but study of sustenmce of reaction needs consideration of spacial 
distribution. If ini tia.lly present poisons were distributed in 
such a way that noar boundaries of reactor their content was much 
higher thrin at centre than reaction would have sustained for 
longer time as well as larger fluence. Such situa.tions coupled 

with continued mineralisation could probo.bly lead to achievement 
21 2 

of even 10 n/cm flu once. 

^<,4 fiecommendat ions for further study :■ 


(1) Kinetic aspects of Oklo reactors should be modelled 

p.nd embeded into criticality and burnup /buildup scheme 
developed in this study. This makes flux as an determinate 
quantity. If this is done we will have a closed loop v/ith 
^eff floix and flux determining tonperature and 

tcaiperature inturn changing For modelling of this 

closed cycle, along with reactor kinetic equations heat 
transfer properties of Okie rock will li£ive to be idealised. 

If all these additions can be done we can simulate Oklo 
reactors on computer. With this other than modifying results 
of sustenance of reaction mode of reaction can be studied. 
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(2) As was pointed out in Sec, 4.2 Oklo moderator has 

interesting fcatu'res. If points to uew concept of safe 
reactor in whic>i neutron leoknge is less and instead it 
could produce useful radioisotopes. It might he an useful 
exerciso to develop this concept of safe reactor. An 
extensive study ccri ho done to find out possible and useful 
modoz-ating material along these lines. The search for 
such moderator will depend on selection of elements which 
crni chemically t orrn dense solid and have high moderating 
power. Choice of these elements can bo such that their 
isotopes can' be used for numerous radioisotope applications. 
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1.2 Thermal data for isotopes of Uraniim and Plutoninmo 
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Contd. . .Appendix- I 

1.4 Radioactive decay constants; 



Pu^^ 



jl35 

Xel35 ; p„149 

i ! 

i ! 

' / 

■ 1 

9.03273x10"^^ 

' -12 
: 3. 37x10 ^ 

1. 67xio"^ i 

( 

2.874x10"^ i 

>. ^ ! 

^ ! 

;2 11x1C‘^;3.6x10" 

1 


1.5 Additional data used s 

Ti25 = 


Value is of infinite resonance absorption integral for U-235* 

(1) Data is taken frcm ML-5800. 

(2) Data taken from Table (18), 2-13, Nuclear Engineering 
Handbook edited by Harold Etherington. 

(3) Miscroscopic transport cross-section for elements other than 
water were calculated by following equation 

0^ = a + (1 - lI^) o 

tv a ^ s 

3 £ 5 € 


All the cross-sections at neutron speed of 2200 m/sec. 
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